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PROBING HIGH-REDSHIFT GALAXY FORMATION AT THE HIGHEST LUMINOSITIES: 
NEW INSIGHTS FROM DEIMOS SPECTROSCOPY 

Kyoung-Soo Lee 1 , Arjun Dey 2 , Michael C. Cooper 3 , Naveen Reddy 4 , Buell T. Jannuzi 5 

ABSTRACT 

We present Keck DEIMOS spectroscopic observations of the most UV-luminous star-forming galax- 
ies at rcdshifts 3.2 < z < 4.6. Our sample, selected in the Bootes field of the NOAO Deep Wide-Field 
Survey, contains galaxies with luminosities of L* < L\jy < 7L* and is one of the largest samples to 
date of the most UV-luminous galaxies at these redshifts. Our spectroscopic data confirm 41 can- 
didates as star-forming galaxies at 3.2 < z < 4.6 and validate the relatively clean selection of the 
photometric candidates with a contamination rate of 11 — 28%. We find that the fraction of Lya 
emitting galaxies increases with decreasing UV luminosity. None of the 12 galaxies with Mtjv < — 22 
(i.e., Ltjv > 3L*) exhibit strong Lya emission. We find strong evidence of large-scale outflows, trans- 
porting the neutral/ionized gas in the interstellar medium away from the galaxy. Galaxies exhibiting 
both interstellar absorption and Lya emission lines show a significant offset between the two features, 
with the relative velocity of 200 — 1140 km s _1 . We find tentative evidence that this measure of 
the outflow velocity increases with UV luminosity and/or stellar mass. The luminosity- and mass- 
dependent outflow strengths suggest that the efficiency of feedback and enrichment of the surrounding 
medium depend on these galaxy parameters. We also stack the individual spectra to construct com- 
posite spectra of the absorption-line-only and Lya-emitting subsets of the UV luminous galaxies at 
z ~ 3.7. The composite spectra are very similar to those of lower- redshift and lower-luminosity LBGs 
samples, but with some subtle differences. Analyses of the composite spectra suggest that the UV 
luminous LBGs at z ~ 3.7 may have a higher covering fraction of absorbing gas, and may be older 
(or have had more prolonged star-formation histories) than their lower- redshift and lower-luminosity 
counterparts. In addition, we have discovered that five galaxies in the sample belong to a massive 
overdensity at z = 3.78. Finally, two galaxies each show two distinct sets of interstellar absorption 
features. The latter may be a sign of a final stage of major merger, or clumpy disk formation. These 
systems are not expected in our sample: their presence implies that frequency of such sources among 
our luminous z ~ 3.7 LBGs may be an order of magnitude higher than in lower redshift and lower 
luminosity samples. 

Subject headings: cosmology:observations - galaxies:distances and redshifts - galaxies: evolution - 
galaxies : formation 



1. INTRODUCTION 

Deep multi-wavelength space- and ground-based 
imaging surveys have enabled the identification and 
study of large numbers of high- redshift galaxi es (e.g., 
iSteidel et al.ll2003L [2001 iGiavalisco et al.ll2004b[ ). These 
studies have made significant strides in character- 
izing various global statistics of the high-redshift 
galaxy population, including their UV luminosity func- 
tion, stell ar mass fun c tion, and clustering proper- 
ties (e.g.. ISteidel et all 119991: IGiavalisco et ail 12004a!: 
Bouwens et alJl2007t IReddy fc Steidelll2009t lOuchi et all 
2004allbl : iLee et all 120061 . I2009L l2012bt iGonzalez et al.1 
20111) . In addition, the overall distribution of dust 
content, stellar population ages, and s izes have been 
determined at different cosmic epo c hs (iFerguson et al l 
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iFinkelstein et alJl2012f >. 

Nevertheless, several questions remain about the origin 
and evolution of these galaxies. For example, what phys- 
ical processes govern the star-formation histories of these 
galaxies? Do the star-formation, assembly and chemical 
evolutionary histories vary as function of galaxy lumi- 
nosity, halo and stellar mass, and environment? If they 
do, what drives these differences? 

While useful constraints have been placed on 
some of these questions ba s ed on statistic al stud- 
ies (e.g., ISteidel et all 12005 ILee et all |20T1 I2012bt 
iPapovich et al.ll20lHlReddv et al .Il2012bl) . the deep spec- 
troscopy necessary to measure the age and metallicity 
of the stellar populations, interstellar abund ances, and 
large - scale outflows , has been lacking (but see lGiavaliscol 
2002; Sh a^levll20TTL for review) . The main challenge has 
been the faintness of most high-redshift galaxies, which 
precludes high signal-to- noise ratio (S/N) spectroscopic 
studies. The usefulness of deep spectroscopy in shedding 
light on the chemical, kinematic, and geometrical details 
within the galaxy has been amply demonstrated in a 
few studies of very rare, bright galaxies (either intrinsi- 
cally l uminous or gravitationally amplified; e.g., MS 1512- 
cB58: iPettini et al.|[2000l [2002T ISteidel et alj|2010h . and 
also by studies of high S/N composite spectra created 
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by averaging many low-S /N spectra (e.g.. IShaplev et al.l 
[2001 ISteidel et alJl200l l201Qt IJones et alJl2012ti . These 
pioneering studies have shown that galactic outflows are 
powerful enough to drive out as much gas as that con- 
sumed by the star-formation at velocities close to the 
escape velocity, t hereby effective l y enri ching the circum- 
galactic medium. IShaplev et al.l ([2003D showed that the 
interstellar medium (ISM, hereafter) in these galaxies 
must be clumpy based on the observed anti-correlation 
between the Lya and interstellar line equivalent widths. 
When a galaxy is observed through optically thinner re- 
gions, more Lya photons escape while less interstellar 
absorption will occur. 

Most of these results pertain to galaxies in a rela- 
tively narrow range of UV luminosity L*) and red- 
shift (typically z ~ 2 — 3). This luminosity and redshift 
ranges represent a "sweet spot" for efficient studies of rel- 
atively large numbers of UV-luminous galaxies, since the 
space density of high-redshift galaxies falls precipitously 
at L > L* while the success rate of spectroscopic identifi- 
cation declines to lower luminosities and higher redshifts. 
However, several studies have begun to suggest that 
galaxies of different luminosities may exhibit different 
characteristics (e.g., the fraction of Lya emi tting galax- 
ies or the stellar- mass- to-UV-light ratio: see lStark et al.l 
120101 iLee et aLll2012bD . It is therefore important to in- 
vestigate the full luminosity range of high-redshift galax- 
ies and determine how representative the current phys- 
ical picture, offered by spectroscopic efforts of L > L* 
galaxies, is to more- or less- luminous galaxies and/or at 
higher redshift. Widening the dynamic range will also 
improve our understanding of the physical processes of 
the feedback and mass assembly histo ries within different 
dark matter ha l o environments (e.g.. | Ouchi et al.l l2004bt 
ILee et alJl2006t iHildebrandt et al.ll2()()7t). 

In this paper, we present a pilot study of spectroscopic 
observations of higher-redshift (3.2 < z < 4.6), higher- 
luminosity (L* < L < 7L*) galaxies with the primary 
aim of critically assessing our current understanding of 
galaxy formation at high redshift. 

We use (0,0^ 8,^100) = (0.27,0.73,0.8,0.7) 
dKomatsu et al.ll2011). Ma gnitudes are given in the AB 
system (IQke fc GunnT[T983l) unless noted otherwise. 



2. DATA AND PHOTOMETRIC SELECTION OF GALAXIES 
AT Z ~ 3.7 

Our photometric candidates are selected from the 
Boote s field of the NOAO D eep Wide-Field Survey (ND- 
WFS; Uannuzi fc Devi[l999t ). Briefly, the NDWFS B W RI 
band data were obtained using the Mosaic camera on the 
Mayall 4m Telescope, and consist of 27 separate point- 
ings covering a contiguous 9.3 deg 2 area. 

For the spectroscopic target selection, we focused on 
three (3) pointings obtained under the best observing 
conditions. We chose pointings where: (1) all three pho- 
tometric bands have seeing better than 1" and within 
0.1" of one another; and (2) the 50% completeness limit 
in the three bands are deeper than Bw > 26.8, R > 25.8, 
and I > 25.1 (Vega)0. The completeness limit was de- 
termined by adding artificial stellar objects (convolved 

6 The NDWFS images are calibrated in Vega magnitudes. The 
correction to be applied to convert to the AB system is 0.019, 
0.215, and 0.459 mag for the Bw, R and /-band, respectively; e.g., 
Iab = -fvega + 0.459. 
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Fig. 1. — The locations of the NDWFS sources on the Bw — 
R-R — / plane are shown together with the Lyman-break color 
selection criteria outlined by the black solid line. The colors are 
based on Vega magnitudes. The spectroscopic targets are overlaid 
in color showing the interlopers (z < 3.2; purple), LBGs at 3.2 < 
z < 4.2 (blue), LBGs at z > 4.2 (green), and the sources that could 
not be identified (red). 



with a Moffat profile matched to the measured seeing) 
to the images and recovering them with SExtractor us- 
ing the identical parameters as the real data (tabulated 
on the survey web page)Q These requirements allow ro- 
bust photometric color selection without degrading the 
images by convolving them to a common seeing, and a 
maximum surface density within the multi-object masks 
down to / pa 25.0. 

Our photometric sample is selected by applying a 
Lyman-bre ak color-selection t e chnique to the BwRI 
data (e.g., ISteidel et al.l 119991: iGiavalisco et al.l I2004at 
IBouwens et al.ll2007f) . The Lyman-break technique is de- 
signed to identify a UV bright star-forming galaxy (or 
Lyman Break Galaxy; LBG, hereafter) with a strong 
spectral break (at a rest-frame wavelength A < 1216 A) 
resulting from absorption by the intervening Lya forest. 
At 3.3 < z < 4.3, this break lies between the Bw and 
R bands. The color criteria (in Vega magnitudes) used 
to select the candidates are identical to that used for 
a larger sample of photometric candidates described in 
ILee et all (|20ll . 



(B w -R) > 1.2 + 3.0 x(R-I) 
(R-I)> -0.3 S/N(R) > 3 



(B w -R)> 2.0 
S/N(I) > 7 (1) 



The total number of sources satisfying the color selection 
criteria are 738. In Figure [TJ wc show the color locations 
of all sources together with the selection window. The 
spectroscopic targets are shown as colored dots; these are 
discussed further in fJ3] 

3. KECK DEIMOS OBSERVATIONS 

Spectroscopic observations of our Lyman break galaxy 
candidates were carried out on the night of May 10 
2010 using the DEep Imaging Multi-Object S pectrograph 
(DEIMOS) on the 10-m Keck II telescope (|Faber et all 
120031) . We used 6000 line/mm grating blazed at 7500 A, 

7 http : //www . noao . edu/noao/noaodeep/DR3/dr3-descr . html 
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Fig. 2. — (Left:) The magnitude distribution of our targets (66), confirmed LBGs (41) at z = 3.2 — 4.7, and unidentified sources (20) 
are shown as a function of /-band magnitude. (Right:) The redshift distribution of the spcctroscopically confirmed LBGs and QSOs in 
the NDWFS Bootes field. Also shown are the redshift distribution of three QSOs and five galaxies in the large-scale structure at z = 3.78. 
The expected redshift distribution from our photometric simulations is shown in dotted line. 

we were unable to determine redshifts, 2 show feature- 
less red continua, indicating that they are likely dusty 
galaxies at z < 2. An additional 5 galaxies show a 
clear continuum break consistent with being a Lyman 
break galaxy, but the signal-to-noise is too poor to iden- 
tify emission or absorption features to confirm (or mea- 
sure precise) redshifts. Finally, no information could be 
obtained on 2 sources because the slitlet was contam- 
inated and therefore masked out during the reduction. 
Hence, the estimated contamination rate from the spec- 
troscopic identifications is 11% (5/46), where we have 
included the 3 < z < 4 QSOs as contaminants. The 
lower limit on the contamination rate for the full sam- 
ple is also 11% (= (5 + 2)/(66 - 2)) where the total of 64 
galaxies excludes the two sources which were masked out. 
In the most pessimistic case where all of the unidentified 
sources are not LBGs, the contamination rate is 28%; 
1 — (41 + 5)/(66 — 2). The magnitude distributions of 
our targets, confirmed LBGs, and unidentified sources 
are presented in Figure [5] (left). 

For the majority of confirmed LBGs (36/41; i.e., 88%), 
a clear UV continuum and at least two interstellar lines 
were detected with sufficient S /N to determine interstel- 
lar redshifts (Table [2]). The remaining 5 LBGs are iden- 
tified only by their Lya emission. Of those 5, we detect 
a UV continuum in three galaxies but at a low S/N that 
we could not measure interstellar redshifts. Therefore, 
the confirmed LBGs at z < 4.2 represent approximately 
a (continuum-)flux-limited sample minimally affected by 
the spectroscopic bias which preferentially identifies line 
emitters. 

Figure [2] (right) shows the redshift distribution of all 
spectroscopically identified sources at z > 3, including 
41 LBGs and 3 high-redshift QSOs. The observed high- 
redshift tail (at z > 4.2) can be explained by the spec- 
troscopic selection effect that favors redshift identifica- 
tion for Lya emitting sources. Of the seven sources 
at z > 4.2, four are secondary targets, detected at 
a low S/N in the /-band, and five are observed with 
a high Lya rest-frame equivalent width (Wo > 18 A). 
The average equivalent width for those five sources is 
Wo — 24.4A, much higher than that measured for the 
full sample (Wo = — 1.2A). Among the measured red- 
shifts are five galaxies lying within a small redshift range 
(z = 3.782 — 3.787) and within 1.5' from each other (all 



yielding a spectral coverage of 4300 — 9600A at the dis- 
persion of 0.65A/pixel. 

We designed 3 masks (BDm_01, BDm_02, and 
BDm_03) covering a total of 66 candidates. The mask 
locations were selected with the aim of maximizing the 
number of Lyman break galaxy candidates within each 
DEIMOS field of view. The photometric candidates were 
put in two categories as "primary" (S/N (I) > 10) and 
"secondary" (7 < S/N(I) < 10) targets, with the pri- 
mary targets driving the DEIMOS field placement. Ta- 
ble Q] summarizes the basic properties of each mask. The 
slit width was 1".2 in all three masks. The effective spec- 
tral resolution estimated from the sky lines is 3.2A. Each 
mask was observed with 2.0 — 2.5 hours in good con- 
ditions, with seeing between 0".6-0".8 throughout the 
night. 

The data were reduced using the spec 2d IDL pipeline 
developed for the D EIMOS instrument (jNewman et al.1 
120121 iCooper et al.1 12012T ) . The wavelength calibration 
was carried out using observations of Ne-r-Ar+Kr+Xe 
arc lamps. The two-dimensional sky-subtracted spec- 
tra were inspected visually to validate the automatic ex- 
traction of the one-dimensional spectra by the DEIMOS 
pipeline. There are a few cases where another galaxy 
fell into the slit and affected the sky subtraction, or 
the galaxy was located too close to the edge of the 
slit. In these cases, re-extraction using the IRAF routine 
noao .twodspec . apall significantly improved the qual- 
ity of the ID spectra. 

Redshifts were determined by visual examination from 
both one- and tw o-dimensional spectra. We used the 
SpecPro software (iMasters fc Capakll2011[ ) to check the 
spectra simultaneously with the available imaging data. 
Redshift identification is made by looking for the pres- 
ence of Lya emission, interstellar absorption lines, and 
the continuum break due to the Lyman alpha forest. Of 
the 66 Bw-band dropout candidates observed during the 
run, we have determined redshifts for 46 sources. In Ta- 
ble[U we list the basic properties of the 46 galaxies/QSOs 
including their coordinates, redshifts, and magnitudes. 

Of the 46 spectroscopic identifications, 41 sources are 
galaxies at the expected redshift range, 2 are low-redshift 
interlopers at z = 0.373 and 0.631 (identified via [On] 
3727 and Ha emission lines), and 3 are QSOs at z = 
3.289, 3.835, 4.041. Of the remaining 20 sources for which 
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TABLE 1 

Summary of DEIMOS Multi-object Slit Masks 



Mask 


Field a 


RA b 


DEC b 


N primary 


-^secondary 


<cx P (hr) 


BDm_01 


NDWFSJ1426p3236 


14:27:10.30 


32:34:11.0 


20 


9 


2.50 


BDm_02 


NDWFSJ1431p3236 


14:31:26.41 


32:25:40.8 


14 





2.25 


BDm_03 


NDWFSJ1437p3347 


14:37:34.37 


33:36:36.1 


16 


7 


1.99 


All 








50 


16 





a The name of the fields in the NDWFS survey webpage 
The J2000 coordinates of the center of each pointing 



TABLE 2 

The Summary of the Spectroscopic Sample 



ID 




RA 






DEC 




I mag 


Class a 




2IS C 


W),Lya d 


UV continuum 


FLAG([3.6 


BD13355 


14: 


,27: 


15. 


.864 


32 


:32 


:34 


.66 


24. 


.27 


P 




0.373 




Y 


N 


BD6909 


14: 


27: 


12. 


.180 


32 


:25: 


:58 


.69 


25 


.08 


P 




0.631 


~ 


Y 


N 


BD28913 


14: 


:27: 


08. 


.131 


32 


:41 


:08 


.92 


23 


.38 


P 


— 


3.289 


[QSO] 


Y 


Y 


BD106241 


11: 


31: 


06. 


.041 


32 


:30: 


:11 


.12 


22 


.24 


P 


— 


3.835 


[QSO] 


Y 


Y 


BD10334 


11: 


27: 


1. 


.982 


32 


:29: 


:30 


.84 


23 


.78 


P 


4.041 




[QSO] 


Y 


Y 


BD97772 g 


14: 


31: 


17. 


057 


32 


:22 


■56 


.21 


25 


.06 


P 


3.782 




5.1 


Y 


N 


BD98176 g 


14: 


31: 


48. 


811 


32 


:23: 


:23. 


.10 


24 


.61 


P 


3.784 


1 


14.2 


Y 


N 


BD99073 


11: 


31: 


37. 


.656 


32 


24: 


:17 


.39 


24 


.16 


P 


3.730 


3.717 


0.9 


Y 


IIC 


BD99847 g 


14: 


31: 


37. 


711 


32 


:25: 


:03. 


.07 


24 


.54 


P 


3.783 




21.3 


Y 


Y 


BD100006 S 


14: 


31: 


5 1 


209 


32 


25; 


:11. 


.06 


23 


.76 


P 


3.787 


■ 


10.9 


Y 


Y 


BD100902 g 


1 1: 


31: 


37. 


001 


32 


:25 


:57. 


.54 


24. 


.45 


P 


3.782 


3.777 


24.9 


Y 


N 


BD102774 


1 1: 


31: 


.15 


.130 


32 


:27. 


:28 


.91 


24. 


.10 


P 




3.753 


-3.8 


Y 


Y 


BD98449 


1 i: 


31: 


.12. 


.230 


32 


:23 


:40 


.67 


24 


.75 


P 


I. : 




26.3 


Y 


N 


BD100141 


11: 


31: 


10. 


.255 


32 


:25. 


:18 


.41 


24. 


.72 


P 




3 'V 


-8.0 


Y 


Y 


BD7645 


14: 


■97. 


I -1 


.009 


32 


.ZD. 


. A X 
.10 


.17 


Zl. 


U.) 


P 




3.685,3.679 h 


-7.2 


Y 


Y 


BD8040 


11: 


97- 


]i 


.386 


32 


■97 


■ns: 

UO 


.03 


9^ 


y 1 


S 


; ■ 


3.666 


15.1 


Y 


N 


BD9126 


14: 


.27: 


15. 


.394 


32 


:28 


:15 


.20 


24 


.59 


S 


4.091 


- 


> 19.8 


N 


N 


BD10059 


11: 


:27: 


01. 


.711 


32 


:29: 


:13 


.56 


24. 


.85 


S 




3.240 


-8.8 


Y 


Y 


BD10252 


14: 


.27: 


01. 


.975 


32 


:29 


:25 


.30 


23 


.80 


p 


4.073 


4.062 


9.6 


Y 


Y 


BD11113 


14: 


.27: 


08. 


.621 


32 


:30 


:21 


.67 


24 


.26 


p 


3.659 


3.653 


5.9 


Y 


N 


BD11352 


11: 


:27: 


.08. 


.940 


32 


:30: 


:35 


.24 


24. 


.20 


p 




3.742 


-5.2 


Y 


Y 


BD15308 


14: 


.27: 


01. 


.745 


32 


:34 


:27 


.66 


23 


.90 


p 




3.301 


-7.6 


Y 


Y 


BD13266 


14: 


:27: 


09. 


.893 


32 


:32 


:28 


.79 


23 


.33 


p 




4.062 


-4.2 


Y 


Y 


BD18449 


14: 


26: 


57. 


.403 


32 


:38 


:19 


.03 


22 


.56 


p 


3.744 


3.726,3.744 h 


-2.6 


Y 


Y 


BD19284 


11: 


:27: 


00. 


.732 


32 


:39: 


:24 


.66 


24. 


.31 


S 




3.260 


-3.7 


Y 


Y 


BD28630 


11: 


.27: 


05. 


933 


32 


:41: 


:23 


.46 


24. 


.98 


S 


3.844 


3.839 




Y 


N 


BD29677 


14: 


26: 


59. 


.578 


32 


:40 


:31 


.76 


24 


.76 


s 


4.205 


4.195 


33.2 


Y 


Y 


BD8859 


14: 


:27: 


05. 


.875 


32 


:27: 


:56 


.74 


23 


.25 


p 




4.099 


-3.8 


Y 


Y 


BD12027 


11: 


26: 


56. 


.448 


32 


:31: 


:U 


.81 


24. 


.74 


p 




4.070 


-6.1 


Y 


Y 


BD12510 


14: 


.27: 


01. 


.061 


32 


:31 


:41 


.15 


24 


.16 


p 


4.717 


4.704 


53.9 


Y 


Y 


BD13763 


11: 


:27: 


13. 


.939 


32 


:32 


:56 


.15 


24. 


.80 


p 


3.672 


3.665 


6.7 


Y 


UC 


BD14778 


11: 


.27: 


06. 


.732 


32 


:33: 


:59 


.69 


24. 


.70 


p 


3.987 


3.980 


-0.5 


Y 


uc 


BD235717 


14: 


37: 


02. 


.976 


33 


:32 


:30 


.37 


24 


.58 


s 


3.826 


3.814 


31.9 


Y 


N 


BD236105 


11: 


37: 


18. 


.703 


33 


:32 


:50 


.39 


24. 


.79 


S 


4.219 




> 18.3 


N 


N 


BD239105 


14: 


37: 


21. 


.862 


33 


:35 


:34 


.62 


23 


.96 


p 




3.900 


-13.4 


Y 


Y 


BD239303 


14: 


:37: 


38. 


.177 


33 


:35 


:43 


.12 


24 


.86 


s 




3.828 


-1.9 


Y 


Y 


BD239473 


11: 


37: 


.52. 


.464 


33 


:35: 


:50 


.75 


24. 


.59 


p 




3.690 


-2.5 


Y 


PC 


BD239555 


14: 


37: 


.54. 


.758 


33 


:35 


:55 


.39 


24 


.75 


p 




3.230 


-1.2 


Y 


Y 


BD240454 


14: 


37: 


,18. 


.238 


33 


:36 


:46 


.37 


24 


.49 


p 




3.663 


-1.1 


Y 


Y 


BD240992 


11: 


37: 


11. 


.093 


33 


:37: 


:19 


.88 


24. 


.61 


p 




3.650 


-5.9 


Y 


N 


BD244504 


14: 


37: 


.58. 


.382 


33 


:41 


:06 


.86 


24 


.29 


p 




3.724 


-1.9 


Y 


N 


BD235538 


14: 


:37: 


11. 


.563 


33 


:32 


:18 


.46 


24 
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a "P" is for primary and "S" is for secondary targets, based on the S/N in the detection band 
k The Lya redshift is determined from the peak of the Lya emission 

c The interstellar redshift is determined from the peak of the interstellar absorption line 

The rest-frame equivalent width of Lya, only measured for LBGs 
c The UV continuum detected in the spectrum? 

Spitzer IRAC [3.6 fim] band detection: "Y" for detected, "N" for undetected, and "UC" for uncertain 
s The galaxies in the z — 3.78 structure (sec jj5.lt 
" The galaxies with multiple interstellar rcdshifts (sec i|5.2t 
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l-band AB Magnitude 

Fig. 3. — The normalized /-band magnitude distribution of 
our sample is compared to that of t he Vanzella+Stark sample 
HVanzella et alj |2009: Stark ct al. 201Cj). The galaxies m our sam- 
ple are roughly 2.1 times more UV-luminous than those in the 
Vanzella+Stark sample. 

in the mask BDm_02). The close physical proximity of 
these five galaxies strongly suggests that they belong to 
a massive structure (see tj5.1[) . In consideration of these 
factors, the observed N(z) is in a reasonable agreement 
with the redshift distribution expected from photomet- 
ric simulations as shown in Figure [5] (dotted line in the 
right panel; see Lee et al. 2011 for the details of the sim- 
ulations), i.e., that our selection criteria yield a galaxy 
population at z ~ 3.7 ± 0.4. 

4. THE PHYSICAL PROPERTIES OF MOST UV-LUMINOUS 

LBGS 

The spectra provide rich details about the physical 
properties of the most UV-luminous galaxies and show 
striking differences in the galaxy properties at differ- 
ent UV luminosities (i.e., star formation rates). Be- 
fore we proceed, it is useful to compare our sample to 
other large spectroscopic samples (gener ally of less UV- 
lumin o us ga laxies) at the same redshift. IVanzella et al.l 
(|2006l 120091) observed « 50 galaxies at z ~ 3.7 with 
the FORS2 sp ectrograph on the Very Large Telescope. 
More recently, iStark et al.l (|2010[) obtained a large sam- 
ple (« 120) of galaxies at z = 3.0 — 4.5 with the 
DEIMOS on the Keck II Telescope. Both samples 
were drawn from the same photometric candidates se- 
lecte d from the Great Obser vatories Origins Deep Sur- 
vey (jGiavalisc o et al.l l2004bO . using the Lyman break 
technique described in §2. In Figure [3l we compare 
the magnitude distribution of our sample to those of 
the combined Vanzella+Stark sample. The difference in 
the median luminosities between the two samples is 0.8 
mag. As the redshift selection functions for these sam- 
ples are nearly identical, the median UV luminosity of 
our sample is roughly 2.1 times higher than that of the 
Vanzella+Stark sample. 

4.1. The Lya Equivalent Width Distribution 

We measure the Lya equivalent width for each galaxy 
directly from the ID spectra. The average continuum 
level c rc( j is determined by taking a weighted average of 
the pixels in the region corresponding to the rest-frame 
1225 — 1255A (and excluding all pixels affected by OH 



sky lines). The line flux F\^ ya is computed by summing 
the flux in excess of the average continuum level c rc d at 
A res t = 1213— 122lA, and the rest-frame equivalent width 
Wq is then estimated as Wq = i*L yQ /[c rc d(l + z)}AX. 
Positive Wq represents emission. Random errors in the 
equivalent widths are estimated from the uncertainties 
in the line flux and continuum level. 

In Figure IH we show the distribution of the rest- frame 
Lya equivalent widths for all LBGs in our sample. Even 
though a half of our sample exhibit some level of Lya in 
emission, the median equivalent width is negative. For 
the 36 "field" galaxies, the median Wq is — 1.2A. On the 
other hand, the 5 "protocluster" galaxies appear to have 
a much higher median EW of 14. 2A. There are only 7/41 
(17%) galaxies that Wo > 20A in our sample. Of those 7 
galaxies, 4 are at z > 4.2 and additional 2 belong to the 
z = 3.78 structure (right panel of FigureH]). If we accept 
the possibility that the majority of z > 4.2 galaxies are 
selected because of their strong Lya emission, and that 
the galaxies in the z — 3.78 structure are not represen- 
tative of the remaining (z < 4.2) field sample, then only 
1/32 "field galaxies" (3%) have W > 20A (BD235717, 
Wo = 32A) . These e stimates are c onsist ent with the mea- 
surements made by IStark et al.l (|2010L see their Figure 
13) and extend these to higher UV luminosities. 

To investigate the frequency of high-EW sources with 
UV luminosity, we compute the luminosity of the sources 
at the rest- frame 1700A ( Figure |4j right). First, we 
adopted the iMadaul ()1995l ) prescription and calculated 
the effective flux attenuation by the intergalactic medium 
(IGM) in the i?-band as a function of redshift, which 
changes from Am = at z < 3.6, to Am = 0.17 at 
z = 4, to Am = 0.39 at z = 4.3. Second, from the IGM- 
corrected (R—I) color, we estimate the UV spectral slope 
/3 assuming a power- law SED (i.e., f\ cx A' 3 ). The ab- 
solute magnitude at A — 1700A, Mi 700, is computed by 
interpolating from the /-band total magnitude. For com- 
parison, we also sho w the locations of the B dropouts in 
the GOODS South (jVanzella et all 120091 ) as open trian- 
gles in the same figure. We computed the UV magnitudes 
at 1700A, by interpolating from the (i — z) color, both of 
which are free of inte rgalactic absorption a t this redshift 
range. We note that IVanzella et al.1 (|2009l ) did not mea- 
sure the equivalent widths for the galaxies with no Lya 
emission, and therefore the triangle points at zero EW 
should be considered as an upper limit. 

From the combined sample, the increasing frequency of 
higher- EW sources (i.e., those with stronger Lya emis- 
sion) towards faint luminosities is evident from these re- 
sults. Interestingly, we find that there are no strong 
Lya emitters (i.e., with Wo > 10A) at M uv < -22.2 
(L\jy > 3L*^ 3 7 ) in the combined sample. Clearly, bet- 
ter statistics are needed to place stronger constraints on 
the Lya-emitting fraction at the brightest end. 

4.2. Large-Scale Outflows 

Interstellar absorption lines that are blue-shifted 
with respect to the galaxy's systemic redshift are 
frequently observed in local and high-redshift star- 
forming galaxies over a wide range of luminosi- 
ties and masses. These observations provide strong 
evidence of galactic-scale outflows that remove a 
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Fig. 4. — Left: The distribution of the rest-frame Lya equivalent widths for the entire spectroscopic sample (top panel) and separated 
into the 5 galaxies in the z = 3.78 structure and the 36 "field" galaxies (bottom panel). The difference in the median values (Wo = 14. 2A 
and — 1.2 A, respectively) between the two samples is significant. Right: The distribution of the Lya EWs as a function of UV luminosity 
suggests a lack of strong Lya emitters at the very bright end. The majority of the high-EW sources are at the high-redshift tail of our 
selection function, i.e., at z > 4.2, which is likely a result of the spectroscopic selection bias. On the other hand, the relatively strong Lyo 
emission from the z = 3.78 structure appears to be real, with 4/5 sources observed with Wo > lOA. When z > 4.2 galaxies and those at 
z = 3.78 are discarded, only 1/29 (i.e., 3%) galaxies have Wo > 2dA. 
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Rupke et alj|2005b lAdelberger et al.ll2q05HVanzella et all 
2009tlWeiner et al.H2009b iSteidel et al.ll2010D . The quan- 



tity of material ejected from a galaxy into the inter- 
galactic medium and the physical process(es) responsible 
are critical to our understanding of "feedback" , which 
is thought to shape a galaxy's star-formation, chemical 
and evolutionary history, as well as the chemical history 
of the intergalactic medium. In particular, the relative 
efficiency of the feedback in galaxies of different luminosi- 
ties /masse s can shed ligh t on their respective asse mbly 
histories. iMartinl ()2005fl and iWeiner et afl (|2009| ) ob- 
served that the outflow velocity increases with the star 
formation rates of the galaxy (local ultra-luminous IR 
galaxies for the former, an d z ~ 1.4 star- form ing galax- 
ies for the latter; also see iRupke et al. 2005). On the 
other hand, several more recent studies of local and 
interm ediate-redshift gal a xies do not find suc h a strong 
trend (jChen et all 120101 : iKornei et all 120121 ). and the 
question regarding the luminosity dependence of feed- 
back remains debated. 

At high redshift, measurements of outflows based on 
the largest sample of star-forming galaxies at z ~ 2 — 3 
do not show any e vidence of a luminosity dependence 
(jSteidel et all 120101 ) and instead find marginal evidence 
for the opposite trend, i.e., that more massive galaxies 
have a lower outflow velocities as traced by the cen- 
troids of interstel lar absorption lines (see Figure 3 of 
ISteidel et al.ll20 10D. The discrepancy between studies at 
different redshifts is hard to understand, but may be in 
part explained by large scatter inherent to the measured 
outflow velocities, which sensitively depends on the view- 
ing angle of the observer as well as the opening angle of 
the winds with respect to the galactic axis. Further- 
more, different line diagnostics studied at different red- 
shifts may trace different kinematic components in the 



ISM, making direct comparison challenging. For exam- 
ple, the Na D doublet (arising in neutral gas) is used 
to trace the outflowing gas in local galaxies, whereas 
intermediate redshift studies use Fe II, Mg II, Mg I, 
which generally have higher ionization potentials. Fur- 
thermore, galaxies at lower redshift can be (and typically 
are) observed at a higher spectral resolution, where the 
systemic and outflo wing components can be more eas- 
ily separated (e.g., iChen et all l2010l ): this approach is 
impractical for faint, high-redshift galaxies. Finally, the 
limited dynamic range in galaxy properties in most high- 
redshift samples studied to date may also play a role in 
erasing any trend of outflow velocity with other galaxy 
properties, especially given the inherently large scatter 
in the measurements or in the intrinsic distributions. In- 
deed, the local correlation is largely driven by the mea- 
surement s of dwarf starbursts and (Ultra-) lum inous IR 
galaxies (|Schwartz fe Martin|[200llMartinll2005l) . In this 
section, we attempt to improve on previous studies of 
outflows by extending the investigation to the more lu- 
minous high-redshift sample of galaxies uncovered by our 
survey. 

The outflow velocity is, by definition, the velocity of 
the ISM material relative to the galaxy's systemic red- 
shift. The latter is usually measured by either stellar 
absorption lines or non-resonant emission lines from ion- 
ized or neutral gas. Unfortunately, most nebular emis- 
sion lines in the rest- frame UV (e.g., Lya and C IV; cov- 
ered by our observed optical spectra) arise from resonant 
transitions and are strongly affected by radiative trans- 
fer effects; the ones that are not resonant (e.g., He II, 
C III]) are generally too weak to have significant detec- 
tions in the existing spectra. In addition, at the redshift 
of our sample the most prominent nebular lines (such as 
Ha, [O m]4959,5007, [O n]3727) are redshifted into the 
near- infrared (some beyond the if-band) at z > 3.8. 

Instead, we measure the relative velocity between the 
Lya emission and that of the interstellar lines as a proxy 
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Fig. 5. — Left: The red histogram shows the distribution of the relative velocity between Lya emission and interstellar absorption lines, 
measured fro m the 18 galax i es tha t clearly show both. The distribution of less UV-luminous z ~ 3.7 galaxies in the GOODS South field 
measured by Vanzclla et al. (2009) is shown in a blue hatched histogram. The median velocity of the galaxies in our sample is larger by 
f» 160 km s~ 1 than th e Vanzella sample. The same quantities measured for a large sample of z ~ 3 galaxies (Shapley ct al. 2003) and 
that of MS 1512-cB58 (Pcttini ct al. 2000) are also indicated. Right: The correlation between the relative velocity inferred from the Lya 
emission and interstellar absorption redshifts, A(«L ya — f a bs)i an d UV luminosity (M1700) is shown together with the Vanzella sample 
(filled squares). The galaxies that are detected (undetected) in the IRAC data are shown in filled circles (triangles). 

for the outflowing medium. iVerhamme et al.l (2006) 
showed that the emergent Lya profile observed at high 
rcdshift can be qualitatively reproduced by a simple 
model in which the energy from supernovae explosions 
drives an expanding shell of uniform velocity V^p. In 
their model, the Lya photons that are back-scattered 
from the receding end of the expanding shell have a 
higher escape fraction by being scattered out of the core 
of the line; this results in redshifting the peak of the Lya 
emission by a relative velocity of Sv — 2V cxp . On the 
other hand, interstellar absorption in the advancing side 
(i.e., the observer's side) of the expanding medium causes 
a trough that is blueshifted by Sv — — T4 xp . In reality, 
there is observational evidence that some of the details of 
the Lya profile is inconsistent with suc h a simple model 
((Quider et al.1 120091 : iKulas et al.l 12012( 1 not to mention 
that there is evidence that ISM is neither uni form nor 
at a constant velocity (e.g.. iSteidel et aT1l2010( l. Never- 
theless, if both Lya emission and interstellar absorption 
are primarily shaped by the outflowing medium, the rel- 
ative velocity between the two should roughly trace the 
outflow velocity to first order. 

We measured the relative velocity AV = V\, ya — Vi$l 
between the interstellar and Lya features for the 17 
galaxies with spectra of sufficient signal-to-noise ratio 
to detect both. The velocity ranges from 235 km s _1 
to 1140 km s _1 as presented in Figure [5j The median 
(mean) value for the sample is 527 (531) km s^ 1 with 
the standard deviation of 220 km s _1 . In Figure [5] (left), 
we show the relative velocity distribution of the 17 galax- 
ies in our sample, along with the measurements for the 
highly magnified z = 2.96 galaxy MS 1512-cB58 and the 
ave rage values for z ~ 3 and ~ 2 star-forming galax- 
ies (iPettini et al.ll2000t IShaplev et alll2003t ISteidel et al.1 
120101 respectively). Also shown is the ve locity distribu- 
tion o f 16 galaxies at z ~ 3.7 measured by I Vanzella et al.l 
(2009). The median velocity measured for the Vanzella 
sample is 370 km s _1 , «160 km s _1 lower than the value 
for our sample, although both have a substantial scatter. 
To test the significance of the difference between the two 
distributions, we used the Kolmogorov-Smirnov test (K- 



S test, hereafter) in IDL kstwo.pro. The test returned 
the probability P = 0.286 that the two measurements 
are drawn randomly from the same distribution. 

Using all 33 galaxies (from our study and the Vanzella 
et al. 2009 sample) with AV measurements, we inves- 
tigate whether the outflow velocity varies with UV lu- 
minosity and stellar mass. Figure [SJd shows the outflow 
measure AV as a function of UV luminosity, Mi 700- The 
two horizontal lines mark the median value for the two 
samples. The data show a lot of scatter, as might be ex- 
pected if orientation geometry, radiative transfer effects, 
or other physical processes play a role in determining the 
average AV" of a galaxy. However, there is a suggestion 
that the "upper envelope" of the outflow velocity is larger 
for more luminous galaxies. We used two statistics to test 
for the presence of a correlation: Kendall's r and Spear- 
man p tests. Both methods measure the test statistics 
(psr or tk) and the probability of a null hypothesis (i.e., 
no correlation; Psr or Pk)- The tests on our sample of 
17 galaxies returned (-0.593, 0.012) and (-0.456,0.011) 
for the Spearman p and Kendall r test, respectively, and 
reject the null hypothesis (no correlation) at the ps 99% 
level. Inclusion of the Vanzella sample results in mea- 
sures of (-0.289,0.103) and (-0.216,0.077), slightly weak- 
ening the strength of this statement, but still suggestive 
of a weak correlation between our outflow velocity mea- 
sure AV and the UV luminosity of a galaxy. Because the 
most luminous galaxy in our sample (BD18449) also has 
the highest velocity offset (1140 km s" 1 ), we repeated the 
same tests excluding the object. The probability of null 
hypothesis is 4% (P = 0.04) for our sample (16 galaxies), 
and 17% for the combined sample (32 galaxies). 

At 3.2 < z < 4.2, the IRAC [3.6pm] band samples 
the rest-frame /-band, and thus can be used as a proxy 
for the stellar mass content. Unfortunately, the existing 
data from the Spi tzer Deep Wide-Field Survey (SDWFS; 
lAshbv et~a l. 2009) are too shallow to robustly detect z ~ 
3.7 galaxies with high S/N (the 3a depth at 3.6Mm is 
23.29 AB mag; see Table 5 of lAshbv et alll2009[ ) . Of the 
17 galaxies with outflow measurements, 6 are detected in 
the [3.6^m] SDWFS imaging, 3 are ambiguous, and the 
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remaining 8 are not detected. In Figure [5b, the IRAC- 
detected and IRAC-undetected sources are marked as red 
circles and triangles, respectively. The Vanzella sample 
is drawn from the GOODS-South survey field which has 
much deeper Spitzer coverag e (the 3a depth at 3.6/im 
of 26.75 AB; see Table 1 of iLee et al.1 1201 2b| ). All 16 
galaxies in the Vanzella sample are securely detected at 
[3.6 nm] with magnitudes ranging over 22.88 — 25.43 AB 
mag with a median of 23.96 AB mag. Given the depth of 
the SDWFS data, only two of the galaxies in the Vanzella 
sample would have been detected had they been in the 
Bootes field. 

Within our sample, there is weak evidence of a cor- 
relation between outflow velocity and IRAC brightness. 
The median (mean) relative velocities are AV = (VL yct — 
V abs ) = 622(655)±226 and412(481)±213kms- 1 for the 
IRAC-detected and IRAC-undetected galaxies, respec- 
tively. Removing the outlying galaxy with the highest ve- 
locity, BD18449, from the IRAC-detected subset changes 
the median to 577(566) ± 99 km s _1 . If we include 
the two galaxies that are marginally detected in IRAC, 
the relative velocity decreases slightly to 577(593) ± 219 
and 526(524) ± 108 km s _1 with and without the source 
BD18449, respectively. In short, IRAC-bright galaxies 
appear to have a higher value of AV than their fainter 
counterparts. Because we do not have robust IRAC pho- 
tometry for the majority of our galaxies, we are unable to 
directly test a bivariate correlation. A K-S test indicates 
that the distribution of the IRAC-undetected galaxies in 
our sample is statistically identical to that of the Vanzella 
sample. On the other hand, the IRAC-detected galax- 
ies are likely drawn from a different distribution, with 
just 7.9% probability of a null hypothesis when compared 
against the Vanzella sample. 

When we consider only the galaxies within the Vanzella 
sample, we find no correlation with luminosity, IRAC 
brightness, or stellar mass. The test statistics are incon- 
clusive, and can only rule out the null hypothesis at the 
~ 45% probability level. 

It is not clear why our sample shows a correlation while 
the Vanzella sample does not. One possible explana- 
tion is that our sample probes a wider dynamic range of 
UV and optical rest-frame luminosities and thus is bet- 
ter suited to discerning a trend. The unknown viewing 
angle of our observations relative to the direction of the 
outflows will result in a large scatter in the observed ve- 
locities of outflows, which can wash away any trend if 
the dynamic range is not large enough. While better 
sampling of galaxies at the most luminous or massive 
end is clearly needed to place the constraints on a firmer 
footing, our results are qualitatively consist ent with stud- 
ies of lower-redsh i ft star-forming gala xies ([Martini 120051 : 
iRupke et al .1120051: 1 Weiner et al.|[2 0091 . where the outflow 
velocity increases with a galaxy's SFR and mass. 

4.3. The Average Physical Properties of the Luminous 
z ~ 3.7 Galaxies 

Creating a high S/N composite spectrum, by aver- 
aging over many individual spectra, can allow us to 
discern spectral features that are not typically visible 
in low-S/N individual spectra, and thereby infer the 
physical properties o f the "typica l " sample galaxy in 
a greater detail (e.g., iSteidel et al.l I200H : iShaplev et al.l 



I2003t I Vanzella et all [20091) . The main challenge of con- 
structing a composite spectrum is the fact that we do 
not have a good estimate of the systemic velocity of each 
galaxy. The spectral features we use for redshift identifi- 
cation, Lya emission and/or interstellar absorption lines, 
are biased with respect to the systemic velocity. As dis- 
cussed in i j4.2[ the interstellar absorption lines are often 
observed to be blue shifted with respect to the galaxy's 
systemic redshift, whereas the peak of the Lya emission 
is redshifted due to radiative transfer effects. 

Deriving the systemic redshift using a statistical cor- 
rection is difficult to justify. Most derived corrections are 
based on Lyman break galaxy sa mples at lower luminos- 
ity and / or lower redshift (e.g.. lAdelberger et alj 12003 
ISteidel et aLll2010D . In light of our findings in jL2]that 
the outflow velocity may depend on the galaxy's lumi- 
nosity, the suggested values from those studies may not 
reasonably apply to our sample. Instead, we compute a 
fiducial redshift for each galaxy and use it to de-redshift 
each spectrum. For the galaxies with only interstellar 
absorption lines, which we refer to as "abs"-type galax- 
ies, the fiducial redshift is computed by applying a con- 
stant offset v to the measured interstellar redshift; i.e., 
Zfld = (— v / c+ zjs) / (1+v/ c) where c is the speed of light. 
We assumed the offset value of v = —250 km s _1 , which 
places S v photosphcric absorption line at the correct 
vacuum wavelength, i.e., A = 1501. 76A. As for the galax- 
ies with both Lya and interstellar redshifts ("comp"- 
type), we take the mean of the two as fiducial redshift. 
This approach should minimize the effect from the range 
of outflow velocities observed in our sample on our com- 
posite spectrum, by making a larger correction if the out- 
flow velocity of a given galaxy is larger. Furthermore, the 
mean correction applied to the "comp" galaxies should 
be v ~ —250 km s" 1 with respect to the interstellar red- 
shifts (see Figure the median velocity offset between 
Lya and IS components is 530 km s _1 ), similar to that 
applied to the "abs" galaxies. 

Because we base our stacking on the interstellar red- 
shifts, we therefore limit our analyses only to the galaxies 
with measured interstellar redshifts. All but four galax- 
ies in our sample have reliable interstellar redshifts de- 
termined from at least two absorption lines, while only 
half (20) have detectable Lya emission. Another advan- 
tage of limiting the analyses to those with interstellar 
redshifts is that the composite spectrum created as a re- 
sult should provide a fair representation of a flux-limited 
sample of high-redshift galaxies with no bias towards the 
Lya line-emitting galaxies. 

We created the composite spectra as follows. First, 
we determine the fiducial redshift for each galaxy as de- 
scribed above. Second, we resample each spectrum to 
a common (near) rest-frame wavelength vector using cu- 
bic spline interpolation (IDL interpol .pro with spline 
option). Third, we normalize each spectrum using the 
average flux density at 1430 — 1470A, where the spec- 
trum is devoid of strong interstellar lines. Fourth, at each 
resampled wavelength element, we average the individ- 
ual spectra weighting by the inverse variance (where the 
inverse variance is derived from the la e rror spectrum 
provi ded by the reduction pipeline; see iCooper et alJ 
2012). For each individual spectrum, the regions con- 
taminated by atmospheric A and B band absorption (at 
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A = 7592 - 7675A and 6886 - 688lA, respectively) are 
masked and excluded from the average. Since the error 
spectrum is background-limited, each galaxy has roughly 
equal weight in the average regardless of its luminosity. 

In Figure [|)J we show the composite spectrum con- 
structed from 36 galaxies in our sample; i.e., all excluding 
five galaxies with no interstellar redshift. In the spec- 
trum, we detect C 111 A1176 line and S v AI502, which 
arise from stellar photospheric absorption. The lines are 
at the correct wavelengths, therefore confirming that our 
final spectrum is in fact at systemic. 

4.3.1. The Galaxies with no Lya emission 

In Figure we show the composite spectrum for our 
z ~ 3.7 sample constructed from 20 sources with in- 
terstellar absorption line redshifts and no Lya emission. 
The spectrum clearly shows not only strong interstellar 
lines such as Si 11 A1260, O i+Si 11 A1303, C 11 A1334, 
Si iv AA1393,1402, Si 11 A1526, C iv AA1548,1550, Fe 11 
A1608, and Al 11 A1670, but also weaker features such as 
Ni 11 A1370, Ly/3, Si III A1206, N v A 1239,1243. 

The top panel of Figure [7] compares our composite 
spectrum with that of the average z ~ 3 LBG from 
IShaplev et all (|2003t) . Specifically, the iShaplev et ail 
(|2003| ) spectrum shown is the composite constructed 
by averaging the 199 galaxies in the bottom quartile 
of the Lya equivalent width distribution (i.e., referred 
to as "Group 1" by Shapley et al.; see their Table 3) 
and exhibits the strongest interstellar absorption. At 
A rC st > 1216A, the z ~ 3 and z ~ 3.7 composite spectra 
closely mirror each other, even down to the small scale 
wiggles which likely result from weak nebular/interstellar 
features. The Lya line morphologies blue-ward of the 
trough are also similar in the two spectra, suggesting 
that the column density and velocity distribution of the 
absorbing gas are comparable. While the interstellar 
line ratios in our sample are generally similar to those 
of z ~ 3 galaxies, the equivalent widths appear to be 
slightly larger in our spectrum (see £14.51 later) . 

Despite the general similarities of our z ~ 3.7 sam- 
ple to the z ~ 3 galaxies, there are also several features 
that may be distinct. Most notably, The C iv P Cygni 
emission, which traces stellar winds from massive stars 
and therefore sensitive to the s tellar initial mass f unc- 
tion and/or population ages (see lPettini et al.ll2000l and 
references therein), may be enhanced for the z ~ 3.7 lu- 
minous LBGs. Both spectra (top panel; Figure [7]) are 
normalized at 1450A and the continuum level at 1550A 
immediately following the C IV emission is not well de- 
termined as can be seen in Figure [7] However, adjusting 
the normalization to a higher value would make the P 
Cygni emission more dramatically different. Although 
the difference is not statistically significant (as can be 
seen from the error spectrum on bottom of Figure UJ), 
the composite spectr a of lower-luminosity galaxies at 
z ~ 3.7 presented in IVanzella et al.l (|2009t third row) 
and IJones et al.l (|2012t fourth row in Figure [7J also ap- 
pear to exhibit similarly strong C iv emission, lending 
support to the possibility that C iv emission may have 
been indeed stronger at higher redshift. 

We also detect N v A1240 P Cygni feature at a lower 
significance. The absorption of the N v doublet is not 
seen in the z ~ 3 Shapley composite, but observed in 



the UV spectrum of MS15 12-cB58 to have W = 0.23A 
(middle panel of Figure 171 iPettini et al.l [20021) . In fact, 
the line morphologies of both C iv and N v in our spec- 
trum are more consistent with that of cB58 than with 
the z ~ 3 composite, even though cB58 is observed with 
a much higher covering fraction (the residual intensity in 
the line cores is nearly zero). Like C iv, N v also traces 
stellar winds from massive O stars and is thus also sensi- 
tive to the galaxy's initial mass function, age, and metal- 
licity. Higher-resolution, higher-S/N data should be able 
to shed light on these parameters. 

We also detect possible weak absorption lines of 
Ni 11 A1370, and Ni 11 A 13 17, which are absent in the 
Shapl ey et al.l (|2003l ) spectrum but are detected in the 
UV spectrum of MS1512-cB58 at a W ~ 0.2A level 
(jPettini et all [2002). Ni is one of the Fe-peak elements 
released by Type la supernovae (SNe) and is found to 
be underabundant in t he cB58 spect r um by a factor of 
3 relative to e.g., Si. iPettini et al.l (|2002f) argue that 
the observed underabundance of the Fe-peak elements in 
cB58 is evidence of a young (~ 300 Myr) population age, 
since metal enrichment by longer-lived, intermediate- 
mass stars has not yet caught up with the elements 
released by short-lived OB stars. The tentative detec- 
tion of Ni in our luminous z ~ 3.7 galaxies may suggest 
that these are "older" (i.e., have larger UV-luminosity 
weighted ages) than the lower- luminosity z ~ 3 galaxies, 
perhaps because the more luminous galaxies have a more 
extended star-forming phase. 

We compare our composite spectrum with the 
IVanzella et alj (|2009f ) "abs" -class composite created from 
21 less-luminous z ~ 3.7 galaxies which lack Lya emis- 
sion (third panel of Figure [7]). The absorption features 
are clearly stronger in our spectrum in most cases. The 
difference appears to be mainly luminosity-dependent 
as both composites have similar Lya absorption EW. 
Our z ~ 3.7 composite also has a redder UV slope 
than the Vanzella composite, consistent with the well- 
know n correlation between luminosity and UV colors 
(e.g., iBouwens etaD 120091 : iFinkelstein et alJl2012f ). As 
in the case of the Shapley composite spectrum, Ni is not 
detected in the Vanzella spectrum. 

Finally, the bottom panel shows th e "abs" -type galax y 
composite from the data presented in I Jones "it al. (2012). 
The median luminosity of 24 galaxies included in the 
spectrum is M\jy ~ —21. Their spectrum has a some- 
what redder slope compared to our composite and the 
Vanzella composite, which may be due to systemic er- 
ror in the flux calibration (T. Jones, priv. comm.). The 
most notable differences arc that their spectrum show a 
higher level of Si* 11 emission (at A = 1265, 1309, 1533A) 
and He 11 emission at 1640A. We discuss this in further 
detail in £14.41 Similar to the Vanzella composite, the 
interstellar absorption lines in the Jones composite are 
generally weaker compared to our spectrum. We com- 
pare the interstellar absorption lines from these spectra 
in further detail in t |4.5l 

4.3.2. Lya-emitting Galaxies vs Non-emitters 

The origin of Lya emission in high-redshift galaxies 
is still not well understood. The observed Lya equiv- 
alent width is a sensitive function of the galaxy's age, 
metallicity, and the amount and relative geometry of gas 
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Fig. 6. — The average spectrum of all 36 confirmed star-forming galaxies at 3.2 < z < 4.6 together with the error spectrum on bottom. 
The spectrum is normalized at 1430— 1470A. Also marked are the vacuum wavelengths of prominent stellar photospheric (dashed), nebular 
(dashed-dot), and interstellar lines (strong lines in solid, weaker lines in dotted lines, respectively). Note that while the expected line 
centroids are labeled, not all these features are detected in the composite spectrum. 
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Fig. 7. — The top panel shows av erage spectrum of gal axies without Ly« emission in our sample (black) compared with the composite 
spectrum of the z ~ 3 galaxies from Shaplcy ct al. (2003) (red). The Shapley et al. spectrum shown h ere is their "Group 1" composite 
which averages spectra of galaxies with Lya EW in the bottom quartile of the overall distribution (sec Shaplcy ct al. 2003, for details). 
The lower panels show the same average spectrum of our sample compared with the spectra of MS15T2-cB58 (second pa nel Pctti nTet al.l 
2002); the z ~ 3.7 composite of lower-luminosity galaxies (third and fourth panel, sec Vanzclla ct al. 2009; Jones ct al. 2012, Respectively). 
Both composites are created from galaxies with no Lyc* emission. 
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Fig. 8. — The composite spectra of the Lya-emitting galaxies (black) vs and non-emitters (green) are shown. Both spectra are normalized 
at 1450A. Interstellar absorption is generally stronger in the non-emitter spectrum, consistent with the expectation of a higher gas column 
density in the ISM for the non-emitters. We also detect Ni II in the non-emitter spectrum bu t not i n the emitter spectrum. Conversely, 
S V is visible in the emitter spectrum, but not in the non-emitter spectrum (see discussions in i|4.3.2l l. 
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and dust in the ISM. In theory, comparing the physi- 
cal properties of galaxies with or without Lya emission 
can provide useful information about which parameter is 
the driving factor of the observed Lya emission in galax- 
ies. Recent s tudies suggest that Lya-emitting galaxies 
are less dusty (iGawiser et al.ll2006b iPentericci et al.ll2007l : 
iKornei et al.H2010D on average than the non-emitters, al- 
though it is possible for ve ry dusty galaxies to ex hibit 
strong Lya emission (e.g., Finkelst ein et al.l 120091 and 
the references therein). Furthermore, different studies 
have conflicting results about how the Lya emission is 
related to the galaxy's p opulation age (jPentericci et al.l 
2009; IKornei et al . 1 12010) . 

One complication that has hindered obtaining a clear 
picture has been that the galaxies observed with Lya 
emission tend to be intrinsically less-luminous in their 
continuum emission. Combined with the spectroscopic 
selection effect (i.e., that it is easier to spectroscopically 
confirm galaxies with strong Lya emission than those 
without), it is challenging to compare fairly the physical 
properties of the emitter and non-emitter populations. 

We directly compare the spectroscopic properties of 
the galaxies observed with and without Lya emission by 
stacking them separately. Once again, we only consider 
the 36 (of the 41) galaxies with the measured interstel- 
lar redshifts. Of these, 17 galaxies exhibit Lya emis- 
sion and 20 do not. One galaxy in the Lya-emitting 
subsample has unusually strong N v emission, and we 
exclude this from the stack. We stacked the remaining 
16 Lya-emitters and 20 non-emitters as described previ- 
ously. The spectrum of non-emitters is identical to that 
shown in Figure [7] 

Figure [8] compares the two composite spectra (normal- 
ized at 1450A). The median (mean) UV luminosities of 
the two samples are comparable: Mn 00 is -21.37 (-21.56) 
and -21.52 (-21.62) for the emitter and non-emitter sub- 
samples, respectively. Similarly, the median (mean) I- 
band magnitudes are also comparable: 24.58 (24.40) and 
24.41 (24.28), respectively. 

It is clear from the spectra that most of the strongest 
interstellar absorption lines, including Si II, O i+Si II, 
C II and Si IVA1394, are stronger in the non-emitter 
composite spectrum. The anti-correlation between the 
equivalent widths of Lya emission and interstellar ab- 
sorption was also o bserved in the z ~ 3 galaxies by 
IShaplev et al.1 (|2003[ ). Shapley et al. argued it as evi- 
dence that the emergent Lya strength and interstellar 
absorption features are jointly determined by random 
sight lines seen through the regions of differing optical 
depths: more optically thick sightlines result in more ab- 
sorption/attennuation of Lya and stronger interstellar 
absorption. 

We also observe several subtle but intriguing features 
that have not been previously reported. The two Ni II 
lines detected in the non-emitter spectrum are absent in 
the emitter spectrum. Although both Ni lines are de- 
tected at a low S/N level, the fact that both appear in 
one spectrum and neither does in the other suggests that 
the trend is in fact real. In addition, S V line is clearly de- 
tected in the emitter spectrum but not in the non-emitter 
spectrum. Ni is one of the Fe-peak elements (Ni, Fe, Mn) 
that are released by Type la SNe, while Type II SNe are 
the main producers of S (together with Si, O, Mg, P). An- 



other Type II element, N v, is visible in the non-emitter 
population only. In the context of standard chemical evo- 
lution models, these observations provide unique insights 
into the recent star formation history of the dominant 
stellar population. Based on these findings, we specu- 
late that the luminous non-emitter population may have 
continued star formation for a considerably longer pe- 
riod of time than the Lya-emitting population, as the 
latter have not had enough time to return metals pro- 
duce d by intermediate-m ass, longer-lived stars into the 
ISM (iPettini et al.l 120021) . S V is no t detected in either 
Pones et all (|2012f ) or TVanzella et al.l (|2009f l emitter spec- 
trum, suggesting that the process responsible for this line 
may also depend on luminosity (or, more likely, age). 

4.4. Emission Line Properties 

We detect several weak Si emission line features in 
the composite spectrum. Si n* emission lines (A = 
1265, 1309, and 1533A) are thought to arise from the fine- 
structure transitions within the outflowi ng gas in the 
inters tellar and circumgalactic medium (Shapley et al. 
2003). Photons absorbed by a ground-state Si n ion is 
re-emitted either at the same energy level or at the cor- 
responding fine-structure transition with roughly equal 
probabilities. In Figure [9j we compare our composite 
spectrum (shown in Figu re [6t with t h ose p resented in 
IShaplev et all (|2003h and Uones et al.l (|2012f l at z ~ 3 
and ~ 4, respectively. The Si n* emission profiles in our 
spectrum are si milar to those in t he Shapley spectrum. 
In contrast, the Uones et alj (|2012t ) spectrum shows clear 
excess Si n* emission, most pronounced at 1533A. Jones 
et al. speculate that the photons scatter resonantly via 
Si II transition throughout the circum-galactic medium 
(CGM) and that the Si II* emission arises from a more 
extended region than the Si II absorption. As a result, 
a fixed width slit might preferentially sample the inner 
parts of galaxies at z ~ 3, but more of the outer re- 
gions at z ~ 4, thus resulting in the observed differences. 
Given the range of galaxy sizes, though, this interpreta- 
tion may be simplistic. 

In the context of this argument, our DEIMOS ob- 
servations use 1.2" slitlets, effe ctively samp l ing 8 .6 kpc 
at z = 3.'ifl. In comparison, Uones et al.l (|2012f ) data 
were taken with 1" slitlets (corresponding to 7.0 kpc at 
z = 3.9). The fact that our data cover a wider area 
and still find that the emission level is comparable to 
that of the Shapley spectrum (i.e., not as s trong as the 
Jones spectrum) suggests that the trend that lJones et al.l 
(|2012f ) found is at least not a redshift-dependent trend. 
Alternatively, it is possible that the galaxies in our sam- 
ple ((Afuv) = —21.4) are considerably larger in physical 
size than those in the Jones sample ((Mqv) = — 21.0), 
and as a result, still only a fraction of emission is en- 
closed within the given aperture. However, galaxy sizes 
are only a mildly inc reasing function of luminosity (e.g., 
iBouwens et al.ll2004h . The detailed information on the 
sizes of galaxies in both samples are required to shed 
more light on the possible explanations of the observed 
differences. 

8 As a reference, the mean half-light radius of z ~ 4 galaxi es mea- 
sured in the rest-frame UV wavelengths is 1.7 kpc IjFerguson et al.1 
12001 ). 
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Fig. 9. — The composite spectrum showing th e re gions around Si II* emission lines and He II is compared with composite spectra of 
z ~ 3 and z ~ 4 LBGs from Shaplcy ct al. (2003) and Jones et al. (2012) respectively. The Si II* emission line profiles are very similar to 
those of the z ~ 3 Shapley et al. (2003) composite, while the [Jones ct al. (2012) spectrum shows more pronounced emission. We do not 
detect He II emission in our composite spectrum, although within margin of error (~ 15% of continuum level), our spectrum is consistent 
with the z ~ 3 Shapley et al. spectrum. 



Another notable emission line in the rest-frame ul- 
traviolet spectru m is He II emission , whic h is clearly 
detect ed in the iShapley et al.1 (|2003| ) and I Jones et"aTI 
(2012) samples. The line profile of the He II emis- 
sion (FWHM ~ 1500 km s" 1 at z ~ 3) suggests that 
the emission has been broadened signific antly by fast, 
dense winds from Wolf-Rayet (W-R) stars (jShaplev et al.l 
2003). W-R stars, as descendants of O stars with masses 
> 20 — 30Af Q , are short-lived, and thus the He II emis- 
sion of W-R origi n is strong only in the initial phase 
of star formation (Schacrcr fc Vaccal 119981 ). Assuming 
constant star formation history and solar metallicity, 
iSchaerer fe Vaccal (fl998l ) estimated that EW(He nl640) 
should drop to zero within « 20 Myr. Assuming subsolar 
metallicity (0.42T©), the He il-bright phase decreases to 
just a few Myr. 

Our data do not show a clear sign of He II emission 
(bottom right of Figure |9j also see Table [3]) even though 
our spectrum is in reasonable agreement with the Shap- 
ley composite within (large) errors. We speculate that 
the lack of strong He II emission may be due to the fact 
that the galaxies in our sample may have more extended 
star formation histories than those in the Shapley and 
Jones sample. While a higher-resolution, higher signal- 
to-noise ratio spectrum is necessary to place more robust 
constraints on the strength of the He n emission in the lu- 
minous galaxies, such an interpretation is qualitatively in 
line with the inferences from other observations discussed 
earlier ( U.3.U £|4.3.2|) that the luminosity -weighted ages 
increase with UV luminosity and/or decrease with Lya 
EW. 



4.5. Dependence of Interstellar Absorption Line 
Strengths on Galaxy Parameters 



The strength of interstellar absorption in the UV spec- 
trum is primarily determined by the kinematics, ionic 
column density, and covering fraction of the absorbing 
gas in the ISM, and therefore provides direct probes of 
the physical conditions therein. The fact that the inter- 
stellar absorption EW correlates so strongly with Lya 
EW suggests the same gas that produces absorption fea- 
tures is also responsible fo r attenuating the Lya emission 
(e.g.. IShanlev et al.l 120031 : IVanzella et al.l [2001 ) . When 
the line is saturated, which may be the case for many 
of the strongest interstellar absorption features at high 
redshift, the equivalent widths primarily depend on the 
velocity dispersion of the gas (measured by deconvolved 
full- width- at-half-maximum) and covering fraction (C/), 
and to a lesser degree, on the ionic column density^ 
Hence, the measurements of such a correlation at dif- 
ferent galaxy properties (e.g., luminosity, mass, color, 
morphology) and at different redshifts can provide valu- 
able insights into the changing physical parameters in 
the ISM among the galaxies. Towards that end, we com- 
pare the measurements of interstellar lines with those in 
the literature. We mainly focus on low-ionization inter- 
stellar lines, Si n A1260, O i+Si n A1303, C n A1334, 
and Si n A1527; and use the average of these four as a 
indicator of interstellar line strength, which we refer to 
as Wlis- These lines represent the strongest interstellar 
lines th at trace cold neutral gas i n the ISM, which were 
used in IShaplev et al.l (|2003D and I Jones etall (|2012[ ). 

We measured the Wlis from five distinct composite 
spectra; [1] the full sample (Figure E]), [2,3] Lya-emitting 

9 The equivalent width of an unsaturated absorption line should 
obey Wis °c Mon! for a saturated line, Wis °c b[\n(Ni oll /b)] ' 5 
while /max / /cont = 1 — Cf where / max is the intrinsic flux density 
at maximum absorption, b is Doppler parameter, and Cf is covering 
fraction of the absorbing gas. 
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TABLE 3 

Equivalent Widths Measured From the Composite Spectrum 



Ion 


Arest a (A) 


W 0Ml (A) 


Wo, emitters (A) 


Wo.non — emitters (A) 


Wo,M uv <-21.4 b (A) 


Wo,M uv >-21.4 b (A) 


H I 


1215.67 


4.5 ±0.2 


15.8 ±0.1 


-4.3 ±0.5 


-4.9 ±0.2 


-3.2 ±0.4 


Si n 


1260.42 


-1.8 ±0.2 


-1.4 ±0.2 


-2.3 ±0.3 


-3.0 ± 1.0 


-2.0 ±0.5 


O i+Si n 


1303.27 


-2.7 ±0.2 


-1.9 ±0.2 


-3.4 ±0.5 


-3.0 ± 1.0 


-3.1 ±0.6 


C n 


1334.53 


-2.3 ±0.2 


-1.6 ±0.2 


-2.7 ±0.3 


-2.0 ±0.4 


-2.2 ±0.4 


Si iv 


1393.76 


-1.4 ±0.1 


-0.9 ±0.1 


-1.8 ±0.2 


-1.9 ±0.4 


-1.9 ±0.3 


Si iv 


1402.77 


-1.0 ±0.1 


-1.2±0.1 


-0.9 ±0.1 


-0.7 ±0.3 


-1.1 ± 0.2 


Si ii 


1526.71 


-1.7 ±0.2 


-0.8 ±0.2 


-2.3 ±0.3 


-2.4 ±0.6 


-2.1 ±0.3 


C iv 


1549.48 


-3.1 ±0.3 


-2.8 ±0.3 


-3.3 ±0.3 


-3.4 ±0.6 


-3.1 ± 0.4 


He ii 


1640.40 


0.3 ±0.3 


-0.2 ±0.5 


0.1 ± 0.5 


0.3 ±0.6 


-1.1 ± 1.0 



The vacuum wavelengths 

The sample only includes the non-emitters in the sample (i.e., Wo, Lya < A) 
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Fig. 10. — Equivalent width of low-ionization absorption lines 
as a function of Lya EW. Large filled circles are from our sam- 
ple divided into two bins as Lya "emitters" and "non-emitters", 
while the large open circle shows the value measured for the full 
sample. Filled upward-poin ting triangles show the measurements 
from the iJones et al.l (2012) z ~ 3.9 samples in three bins accord- 
ing to Lya EW, while downward-pointi ng triangles indicate similar 
measurements at z ~ 3 as published in Shaplcy ct al. (2003). Also 
shown are two luminosity bins (M < —21.4 and M > —21.4) of 
non-emitters in filled diamonds. 

galaxies and non-emitters (Figure [8]), and [4,5] two non- 
emitter samples binned by UV luminosity (divided at 
the median luminosity Muv = —21.4). We measured 
the line EWs directly from the composite spectra; we 
also independently measured the same lines from the 
composite spectra o f the four subsamples presented in 
Shapl ev" 3t al.1 (|2003l ). and confirmed that our measure- 
ments return very similar values (within O.lA) to their 
published values. These results are shown in Figure 
[TOl and tabulated in Table El In the figure, the value 
for our full sample is shown in large open circle, while 
the two subsamples binned by Lya EW are shown in 
filled circles. The measurements for two luminosity (non- 
emi tters) bins are s hown in file d diamonds. We al so show 
the IShaplev et all (|2003l ) and IJones et al.l (|2012f ) points 
as downward and upward triangles, respectively. 

It is clear that Lya-emitting galaxies have much weaker 
interstellar absorption (by more than lA), thus confirm- 
ing the earlier studies. However, it is also true that our 
points for non-emitters and for the full sample lie sys- 
temically lower (i.e., stronger absorption) by ~ 2a than 
the other samples when Lya EW is fixed. As a guide, 
we mark the la range by interpolating the three Shapley 



subsamples (except for the highest W^^ ya bin; Group 
4). All three Jones et al. (2012) points are consis- 
tent with this relation within error. The fact that the 
measurements at z ~ 3 and z ~ 4 are very similar at 
comparable luminositie s ((Muv) = —21: IShaplev et all 
l2003t IJones et"aT1l2012fl rules out that we are seeing red- 
shift evolution of the absorbing gas. Rather, it may be 
that the observed discrepancy is a luminosity-dependent 
effect because the median luminosity of our sample is 
—21.4; i.e., 45% more luminous). We attempted to test 
the luminosity-dependence by splitting the non-emitter 
sample into two bins according to UV luminosity (dia- 
monds in Figure [10]), but were unable to measure any 
significant difference due to the l arge uncerta i nties (each 
composite had only 10 galaxies). IJones et al.l (2012) also 
observed a luminosity-dependent trend in their measure- 
ments; that the low-ionization inter-stellar (LIS) absorp- 
tion was weaker for lower-luminosity galaxies when Wh ya 
is fixed. Similarly, they reported some level of variations 
in LIS absorption on other galaxy parameters such as 
UV spectral slope, half-light radius, UV luminosity, and 
stellar mass (see their Figure 9). The trend that they 
observed is in the same direction as our observations; 
that more-luminous/redder/larger galaxies lie below the 
fiducial relation while less-luminous/bluer/smaller coun- 
terparts lie above the same relation. 

Despite relatively large uncertainties in the Jones mea- 
surements and ours, these results show a hint of chang- 
ing physical conditions within the ISM that depend on 
galaxy properties. Since many demographic properties 
of galaxies correlate well with one another (luminosity, 
color, size, stellar mass, total mass), it is not possible 
to pinpoint what may be the driving factor in producing 
the observed trends. Nevertheless, it is still worthwhile 
to speculate about possible physical scenarios that may 
give rise to the subtle differences observed. As mentioned 
above, our sample contains more lumino us and therefore 
presu mably more massive galaxies (e.g. JLee et al.ll2006l 
l2(TTl , We can therefore expect these galaxies to have 
higher velocity dispersions on average, thus contribut- 
ing to increased equivalent widths for the interstellar ab- 
sorption lines. In addition, the covering fraction of the 
absorbing gas may be higher in more luminous galax- 
ies. Such a scenario would not only imply deeper (more 
saturated) absorption lines in more luminous galaxies, 
but would also explain the observed decrease in the frac- 
tion of Lya-emitting galaxies towards higher luminosi- 
ties (as discussed in ij4.ll and Figure [4j . There is also 
circumstantial evidence that more luminous galaxies are 
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Fig. 11. — Left: The relative positions in angular and redshift 
space of 11 galaxies observed on the BDm_02 mask are shown in 
filled circles. The angular distance is computed from the center 
of the five galaxies; RA=14 h 31 m 44.88 s , DEC= 32° 24' 30.24" 
(J2000). In the inset, we show in open histogram the redshift dis- 
tribution determined from the current spectroscopic sample com- 
piled to date. The yellow shade illustrates the normal distribution 
with <t z = 0.4 centered at z = 3.7. The redshift distribution of the 
galaxies on the BDm_02 is shown in filled histogram. Right: The 
zoom-in on the over density region. The angular positions are now 
shown in units of Mpc (physical). 

more reddened, suggesting a higher covering fraction of 
dust than in less luminous galaxies. If the dust and gas 
distributions are similar, this also supports an increased 
covering fraction of gas in more luminous galaxies. Fi- 
nally, it may be possible that the more luminous galaxies 
have ISM that are more chemically enriched. However, 
this last speculation is difficult to confirm without much 
higher quality (and higher resolution) data. 

5. SPACE ODDITIES 
5.1. A Large Scale Structure at z — 3.78 

The spectroscopic redshifts reveal the presence of a 
candidate large-scale structure at z = 3.78. Five of the 
eleven galaxies on the mask BDm_02 are identified at 
redshifts 3.782 < z hya < 3.787, within 1 Mpc (physi- 
cal) from one another. In Figure QTJ we show the rela- 
tive positions of the eleven galaxies observed in the mask 
BDm_02 in redshift and angular distances, and also com- 
pare with the redshift distribution determined from the 
current compila tion of spectros copy of our photometric 
candidates fsee iLee et al.l [20lTl for further detail). As 
evident from the histogram, the redshift distribution on 
the mask BDm_02 is highly unusual in comparison with 
that of the overall population in our sample. 

We determined interstellar redshifts for four of the five 
galaxies, with the values ranging over Zjs = 3.773 — 3.779. 
The difference between the interstellar and Lya redshifts 
for these sources suggest the relative velocity of 310 — 530 
km s _1 . These values are comparable to the relative 
velocities determined for "field" galaxies, which we dis- 
cussed in gjj Two galaxies (BD100006, BD100902) are 
individually detected in the SDWFS IRAC [3.6 fxm] im- 
age, suggesting that they have larger stellar masses than 
the other three, assuming a similar star formation his- 
tory. 



All five galaxies show pronounced Lya emission as can 
be seen in the two-dimensional spectra (Figure [12]) . The 
mean and median equivalent widths are W = 15.3A 
and 14.2A, respectively, much higher than that deter- 
mined for the full sa mple Wp = —1.2 A. Similar observa- 
tions were made by IKuiper et al.l (|2012D for a z = 3.13 
structure near the radio galaxy MRC 0316-257. Three 
LBGs associated with the radio galaxy at z — 3.13 have 
Lya EWs of W = 14.9, 39.3, and 17.8A. Based on 
the effective EW m easured from the stacked spectrum, 
IKuiper et all (|2012D estimated W = 26 A ± 3.8A for the 
3 protocluster galaxies compared with 7.4 ± 2.lA for 12 
field galaxies. The observed high EW is unlikely to be in- 
fluenced by the spectroscopic selection effect, as all three 
galaxies in the MRC 0316-257 overdensity and four (out 
of 5) of the ones in our overdensity would have been iden- 
tified via absorption line features even if the Lya emission 
was not present. Hence, these observations suggest that 
LBGs in the most massive potential wells may have, on 
average, intrinsically stronger Lya emission than their 
counterparts in the field. 

If the stronger Lya emission observed for these galaxies 
is interpreted as higher star formation rates, our results 
would imply that the galaxies in dense environments 
have, on average, higher SFRs than those in the field; 
the opposite of the SFR-density relation observed locally. 
The "reversal" of the SFR-density relation was alread y 
observed at z ~ 1 (|Elbaz et al.ll2007tlCooper et al.ll200a) . 
where a large number of bright star-forming galaxies ex- 
ist in group environments. Considering the fact that 
the SFR-density relation is likely produced as a result 
of quenching that preferentially takes place in massive 
systems, the reversal of such a relation is not surprising 
at high redshift where most galaxies are still well be- 
low the quenching "threshold" in masses. Accurate mea- 
surements of their stellar masses and extinction proper- 
ties (UV spectral slope as a proxy for reddening in UV, 
or direct measurements in the far-infrared) will provide 
the estimates of SFR and specific SFR thereby shedding 
more light on the nature of their Lya emission. 

Interestingly, the only two galaxies in our sample that 
exhibit doubly-peaked Lya emission both belong to the 
z = 3.78 structure (FigurcDD BD100902 and BD98176). 
While the frequency of s imilar sources at z ~ 2 — 3 is only 
slightly lower 20 - 33% (|Kulas et al.ll2012D than that ob- 
served in our sample, we speculate that it is possible that 
these galaxies may have intrinsically different physical 
parameters that affect the transmission of Lya photons 
(e.g., Doppler parameter, H 1 column density, the veloc- 
ity distribution of the ISM). 

Direct compariso n with other known protoclusters at 
high redshift (e.g [Steidel et al i ;2000: O uchi et alJl2005l : 
IDaddi et al.l 120091 : ICapak et al.l 120111 ) is not trivial due 
to the varying depths and search methods of different 
surveys. The depth of the NDWFS data on which our 
selection of candidates is based allows us to probe down 
to L m L* at this redshift range, where the galaxy sur- 
face density is significantly lower (< 0.1 arcmin -2 ) than 
that of typical galaxies identified as protocluster mem- 
bers in other studies. For example, the 'GN20' proto- 
cluster members in the GOODS South field are identified 
down to L w O IL* in their UV continuum level at the 
same redshift (|Daddi et al.l 120091 ). The surface density 
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Fig. 12. — The two-dimensional DEIMOS spectra of the five galaxies at z = 3.78 — 3.79 are shown. Each spectrum is shifted to align 
the spectra in the observ ed w avelengths. The presence of Lya emission is evident in all cases. The close proximity in both angular and 
redshift space (see Figure \TT\ strongly suggest that these galaxies are physically associated. BD100902 and BD98176 have doubly-peaked 
Lya emission. 

of z ~ 3.7 photometric candidates in the GOODS South 
field is w 9 - 10 arcmkr 2 fsee lLee et ail 1201 2bD . 

The current spectroscopy only covers an 5' xl6' of the 
area in the region, and we are unable to determine the 
physical extent of the structure. Deeper imaging and 
more spectroscopy over a wider area is crucially needed 
to better quantify the uniqueness of the structure and to 
directly compare the physical properties and evolution- 
ary state of is constituents with those of the field galaxies 
(K.-S. Lee et al., in prep). 



5.2. Galaxies with Multiple Interstellar Features 

Among the 41 spectroscopically confirmed LBGs at 
3.2 < z < 4.6, we have identified two sources, BD7645 
and BD18449, that clearly exhibit two sets of interstellar 
absorption lines. The one-dimensional spectra of these 
galaxies are shown in Figure [T3] where prominent inter- 
stellar absorption lines are marked at two different red- 
shifts (cyan and red lines). 

BD7645 is continuum-dominated with Lya in absorp- 
tion. The absorption features are very deep with nearly 
zero residual intensity in Si II A1260, C n A1334, Si iv 
A1393, and C iv AA1548,1550 lines. The redshifts are 
3.685 and 3.679, implying the relative velocity of 380 
km s , the v alue comparable to that reported for z ~ 
3-4 galaxies (IShaplev et alj|2003 IVanzella et all 120091 
see also M4.2|) . 

BD 18449, exhibits a weak Lya emission in combina- 
tion with deep absorption trough. We measure the net 
equivalent width to be Wo, Lya = — 2.2A. The residual 
intensity is «0.5, a half of that observed for the BD7645. 
The C II A1334 and Si II A1526 lines are at z = 3.744, at 
rest with Lya, while O I+Si n A1303 and possibly with 
C iv AA1548,1550 lines indicate an interstellar redshift 
of Zis — 3.726. Interestingly, there is a hint of N iv] 
A1486.5, which is sometimes seen in emission but in this 
case, in absorption. There may be another component, if 
interpreted as a third O i+Si n A1303 blend, at an even 
higher velocity even though we cannot rule out the pos- 
sibility of a chance absorber. The relative velocity of the 
two components in BD18449 implied by the two secure 
sets of interstellar redshifts is 1140 km s _1 , the highest 
observed in our sample and twice larger than the median 
value 530 ± 220 km s" 1 (see AO]! . 

Galaxies with multiple interstellar features are ex- 
tremely rare with only a handful of reported cases; 
Among w 3000 galaxies at z = 1.5-3.5 (LBGs, BX/BM: 



ISteidel et al1 [2003. 2004]), there are 10 such sources iden- 
tified (A. Shapley, priv. comm.). T he most well-known 
source is Q0000-D6 (D6, hereafter; ISteidel et all [20Tol 
e.g., see their Figure 9), the most UV- luminous galaxy 
(1Z = 22.88) in their samples with the exception of 
MS 1512-cB58, which is highly magnified. Similar to 
our sources, D6 shows two distinct sets of interstel- 
lar absorption lines at zis = 2.9561,2.9635, where the 
latter appears to b e at the galaxy's systemic redshift 
(jSteidel et al.l 120101) . D6 has a prominent Lya emis- 
sion which is redshifted with respect to the interstellar 
lines at ZL ya = 2.9692. Another luminous source, Q1623- 
BX453 (1Z = 23.38), shares similar characteri stics to D6 
(z svs = 2.1820, zis = 2.1724, z Lya = 2.1838: lLaw et al.l 
12001 . 

The simplest explanation of this phenomenon is that 
two galaxies in close proximity are mistakenly observed 
as a single source. In this scenario, the continuum lu- 
minosities of the two sources should be comparable in 
order to produce two visible sets of absorption line fea- 
tures. The NDWFS I and i?-band imaging data do not 
reveal multiple components for the two cases here; the 
surface brightness profile is consistent with a single iso- 
lated source in both cases. Thus, if the spectra are pro- 
duced by two distinct galaxies, they must lie within 0.5" 
(i.e., 3.5 kpc in physical scale) of each other, perhaps in 
the last stage of merger. While we cannot rule out the 
possibility of witnessing two final-stage mergers, an al- 
ternate scenario would be that the spectrum is produced 
by two massive rotating clump s perhaps in the pr ocess of 
disk formation as proposed bv lDekel et al.l (|2009ft . High- 
resolution imaging data will shed more light on the na- 
ture of these sources. 

6. DISCUSSION AND SUMMARY 

In this paper, we have investigated the physical prop- 
erties of the most UV-luminous (L > L*) star-forming 
galaxies at z ~ 3.7 based on the rest- frame UV spectra 
of 41 galaxies obtained using the DEIMOS spectrograph 
at the W. M. Keck Observatory. The galaxies in our 
sample populate the exponential tail of the galaxy UV 
luminosity function, and have been under-represented in 
most previous studies due to their low surface density. 
They likely reside in more massive dark matter halos 
and are typically found to have larger stellar masses (> 
a few xl0*°M(T) ) than their less-luminous counterparts 
(|Lee et aljfeoill ). The relatively high S/N spectra af- 
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Fig. 13. — The spectra of the galaxies exhibiting multiple sets of interstellar absorption lines are shown. Both error spectra are in 
units of counts. The locations of prominent absorption lines are indicated by vertical lines at two different interstellar redshifts. The 
interstellar redshifts derived for BD7645 (bottom) are z = 3.685,3.679, with the velocity offset 380 km s _1 . On the other hand, BD18449 
have Zxs = 3.744,3.726, wher e the former is at rest with the Lya redshift. The velocity difference is 1140 km s — 1 , more than twice larger 
than the median value ( 84.2^ . 



forded by their optical brightness have allowed us to 
probe a variety of physical properties of individual galax- 
ies as well as those of the population as a whole. Fur- 
thermore, we have discovered interesting classes of ob- 
jects that are rarely found in lower-luminosity galaxies. 
These discoveries raise interesting questions about how 
similar the physical processes of galaxy formation are at 
different scales (i.e., environments, masses, luminosities). 
In what follows, we summarize our findings and discuss 
their physical implications in further detail. 

The spectroscopy presented here validate the clean and 
efficient selection of photometric candidates, which we 
have identified over a much wider area, 5.3 deg 2 of the 
Bootes field, a nd have used to study the general proper- 
ties of LBGs (|Lee et al.ll201lL l2012a[) . We estimate the 
contamination rate to be «11%, mainly due to lower- 
rcdshift (z < 1) interlopers and high-redshift QSOs (iJH 
Figure^]) ■ Even in the most pessimistic case where we as- 
sume all the galaxies that could not be identified spectro- 
scopically are not at high redshift, the contamination rate 
is 28%. The median luminosity of our spectroscopic sam- 
ple is 2.1 times higher t han previous samples of galaxies 
at co mparable redshifts (Vanzell a et al.l20 09; Star k et al.1 
[20101 see Figure [3]). 

We find that the Lya equivalent width, i.e., the 
strength of Lya line emission relative to the UV contin- 
uum emission, in the galaxies in our sample is generally 
much weaker than tha t found for lower-luminosity galax- 
ies in the literature (jVanzella et al.l 120091 : IStark et all 
2010). The median equivalent width of our sample is 
Wq = — 1.2 A, when we exclude galaxies that are likely 
not representative of our sample; these include galaxies 
at z > 4.2 identified mainly via their strong Lya emis- 
sion despite their faint UV continuum, and five galaxies 
that belong to a large-scale structure at z = 3.78. We 
find only one galaxy with Wo > 20A among the 41 galax- 
ies at 3.2 < z < 4.2. Furthermore, at the brightest end 
(M1700 < —22.2), there appears to be no galaxy with 
a positive Lya EW. The trend of a decreasing number 
of Lya emitting galaxies towards a higher luminosity is 
clear in our data, supporting a changing physical pa- 
rameter that affects the escape fraction of Lya photons 
dramatically. It is important to quantify this dependence 
by increasing the number of spectroscopically observed 
luminous (> 3L*) galaxies. 



The most straightforward explanation of this phe- 
nomenon is that the covering fraction of neutral gas in 
the interstellar medium increases with UV luminosity. 
Such a picture is consistent with both the range of Lya 
EW observed at a given luminosity, and the gradual de- 
cline in the fraction of Lya-emitting galaxies towards 
higher luminosities. O ur recent Hersche l observations 
also support this view ([Lee et al.l l2012al) : by compar- 
ing the UV-derived extinction estimate (inferred by the 
UV continuum slope fx oc A^) to a more direct mea- 
sure of dust extinction (measured by the infrared-to-UV 
luminosity ratio), we find that the UV-based estimate 
under-predicts the true extinction by a factor of 3 for 
the most luminous galaxies (i.e., M1700 < —22.2). The 
two estimates can be brought into better agreement if the 
UV l uminous galaxies foll ow an SMC-like dust reddening 
law (jBouchet et al.lll985l ) rather than the Calzetti law; 
the latter is often used as the sta ndard assumption fo r 
high-redshift star- forming galaxies (Calz etti et al.ll2000n . 
The shallower slope for the extinction with wavelengths 
for the Calzetti law is generally interpreted as resulting 
from a clumpy ISM in the local starburst galaxies, which 
effectively allows more "blue" photons to escape with- 
out attenuation, compared to a more uniform coverage 
of the SMC-like ISM. Although the relative geometry of 
dust and neutral gas is unconstrained for high redshift 
galaxies, there is some evidence that they at least par- 
tially trace each other ( Sha pley et al.ll2003h . If the spa- 
tial distributions of dust and gas are largely decoupled, 
our interpretation would be invalid. We caution readers 
that our Herschel results are based on the photometric 
candidates, the majority of which are not spectroscopi- 
cally verified. Even though the spectroscopy presented 
here seems to confirm the robustness of our sample with 
at most 28% of contamination, we cannot rule out the 
possibility that the results may be largely influenced by 
lower-redshift interlopers. 

Another possibility is that the metallicity of the ISM 
increases with UV luminosity. The fact that interstel- 
lar absorption is generally stronger in higher-luminosity 
galaxies is qualitatively in agreement with this scenario. 
However, the line equivalent widths in our observations 
are a combined product of the column density and ve- 
locity dispersion of the absorbing gas. Without knowl- 
edge of the intrinsic line profiles, it is difficult to disen- 
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tangle these two factors and convert the measured EWs 
to the column density (i.e., relative abundances in the 
ISM). Locally, there is evidence that the velocity dis- 
persion is large r in more UV-luminou s and/or more red- 
dened galaxies (jHeckman et al.l ll99&). which is expected 
if more luminous galaxies have larger dynamical masses. 
Measurements of elemental abundance ratios based on 
higher-S/N spectra may be able to constrain the relative 
importance of these two factors in the future. 

We investigated the relative velocity of the outflowing 
interstellar gas in our luminous LBG sample by mea- 
suring the difference between the redshifts of the inter- 
stellar absorption and Lya emission lines, i.e., AV^ u t = 
Vhya — Va,hs- We find that there is a correlation be- 
tween the UV luminosity and AV^ u t in the sense that 
a higher-luminosity galaxy has a higher observed "max- 
imum velocity" (see £ 14.21 Figure [5|). This implies that 
more UV-luminous galaxies have more powerful outflows 
which can potentially remove cold, neutral gas from the 
galaxy's potential well more effectively. Such a mech- 
anism would have a profound effect on the subsequent 
formation and evolution of the galaxy. On the other 
hand, the gas in less powerful outflows may eventually 
trickle back into the galaxy, thereby only delaying the 
star formation rather than preventing it entirely. 

The observed correlation in our st udy is qualita - 
tively in agreement with that found by Martin (2005); 
iWeiner et al.1 (2009) for local and intermediate-redshift 
star-forming galaxies. In c ontrast, obser v ations of high- 
redshift LBG samples by [Steidcl et al. (2010) suggest 
that the outflow velocity (defined as Kbs) is lower for 
galaxies with larger dynamical masses. It is possible 
that the large scatter in the observed velocity produced 
by random viewing angles has eroded any luminosity- 
dependent trend of the outflow velocity. The trend is 
subtle (e.g., Weiner et al. (2009) found V out oc SFR 03 ) 
and the wider luminosity range probed by our current 
study may have helped identify the trend. We also find 
evidence that the outflow velocity is larger for the IRAC- 
detected (and therefore presumably more massive) galax- 
ies even though the low sensitivity of the current IRAC 
data could not provide reliable estimates of the masses 
for individual galaxies. 

In order to study the average spectroscopic properties 
of the galaxies, we constructed a composite spectrum. 
The composite spectrum of galaxies exhibiting Lya in 
absorption is compa red to similar galaxies at z ~ 3 
(|Shaplev et al.1 l200l gXU Figure [7j. The Lya line 
morphology suggests that the column density and ve- 
locity distribution of the two samples are comparable. 
We find tentative evidence of enhanced emission in ClV 
and N v, both indicative of stellar winds, in our spec- 
trum compared to the z ~ 3 spectrum. If confirmed, 
the results would have interesting implications for the 
IMF and metallicity of galaxies at different luminosities. 
We also marginally detect Ni II, which is absent in the 
z ~ 3 composite. In terms of the wind features and the 
relative strength of the Fe-peak element such as Ni, our 
spectrum may be more similar to that of MS1512-cB58 
than the z ~ 3 spectrum. A possible interpretation of 
the differences is that more luminous galaxies (as rep- 
resented by the composite spectrum of our sample) may 
have been forming stars longer and therefore has a higher 



relative enrichment in the Fe-peak elements (produced 
by intermediate-mass, longer-lived stars) than their less- 
luminous counterparts (the z ~ 3 Shapley et al. compos- 
ite). 

We also compared the composite spectrum of galaxies 
with and without Lya emission at comparable UV lu- 
minosities ( M4.3.21 Figure H]). The interstellar line EWs 
are clearly stronger for the non-emitters, consistent with 
other studies in the literature. The trend likely re- 
flects the combined effect from differing covering frac- 
tions and gas column densities of the two populations. 
However, we have found some evidence that there may 
also be intrinsic differences between the two populations. 
Ni II is detected in the non-emitter spectrum only, while 
S v is seen in then emitter spectrum only. The results 
may imply that Lya-emitting galaxies generally have a 
'younger' luminosity-weighted age than the non-emitters. 
The emerging picture is that the galaxy's age depends on 
both UV luminosity and Lya EW. One simple explana- 
tion is that the covering fraction of neutral gas in the 
ISM increases with UV luminosity, and more luminous 
galaxies tend to be older. 

Two independent observations provide circumstantial 
evidence to this scenario. First, the non-detection of 
He II emission in our composite, compared to its clear de- 
tection in th e composite spectra of less- luminous z ~ 3, 
~ 4 galaxies (jShaplev et alj|2003t Pones et alj|2012f >. im- 
plies that Wolf-Rayet stars, which, as descendants of 
massive O stars, dominate the He n emission at the ini- 
tial period (< 20 Myr) of star formation, are no longer 
present (i j4.41 Figure [§]). Second, the strengths of low- 
ionization interstellar absorption lines are clearly higher 
for more UV-luminous galaxies when the Lya equiva- 
lent w idth is fixed ( §4.5[ Figure [TOJ also see lJones et al.l 
2012). Such effects are likely contributed by a higher 
covering fraction, which determines the residual inten- 
sity of the spectrum, for UV-luminous galaxies although 
we cannot rule out the possibility that a higher column 
density and/or different kinematics of the gas are also 
responsible. 

We have also discovered a possible large-scale struc- 
ture at z = 3.78 where five galaxies in our sample lie 
within 1 Mpc (physical) from one another ( £15.11 Figure 
Hip . These galaxies appear to have enhanced Lya emis- 
sion, (Wo, Lya) = 14.2A, Compared tO (W ,Lya) = -1.2A 

for the field galaxies at the same redshift range (Figure 
fl2land l4l). A similar trend was observed bv lKuiper et al.1 
(2012), but both observations are limited by the small- 
number statistics. We are in the process of obtaining a 
better census of the extent of this structure by survey- 
ing a much wider flanking region around these galaxies 
(K.-S. Lee, in prep). More physical insights require more 
comprehensive, deeper coverage of the region, which is 
currently lacking. Within the hierarchical framework 
of galaxy formation, we argue that the discovery of a 
massive structure at high redshift is not at all surpris- 
ing considering the unprecedented cosmic volume sam- 
pled by our candidates. Other wide-area surveys such 
as CFHTLS, SDSS Stripe 82, Dark Energy Survey (and, 
in the distant future, LSST) should be able to identify 
similar systems in large number, providing a clear view 
of how galaxy formation proceeded in the densest, most 
massive environments. 
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Of 41 galaxies in our sample, two show evidence of 
multiple interstellar components that are widely sepa- 
rated in the velocity space f £)5.2[ Figure [T3|) ; the relative 
velocities between the two components are « 400 km s _1 
for one case, and w 1100 km s -1 for the other. Galax- 
ies with multiple interstellar components are extremely 
rare (< 0.3%) in the existing studies at lower redshift 
(z ~ 2 — 3). Both galaxies appear to have a single com- 
ponent with a smooth surface brightness profile in the 
ground-based optical images, suggesting that if the in- 
terstellar features were produced by two distinct galax- 
ies, they must be within 3 kpc from each other, perhaps 
in the final stages of merging. While we cannot com- 
pletely rule out the possibility of a chance observation 
of two such events based on the small-number statistics, 
we also consider an alternate physical mechanism: that 
the interstellar features originate from giant clumps pro- 
duced by the gravitational instability (e.g., iDekel et al.l 
2009). Regardless of their physical nature, the appar- 
ent high incidence rate of such systems in our sample of 
L > L* galaxies warrants further investigation to obtain 
better statistics. 



The data presented herein were obtained at the W.M. 



Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the 
University of California and the National Aeronautics 
and Space Administration. The Observatory was made 
possible by the generous financial support of the W.M. 
Keck Foundation. We thank the staffs of the W. M. Keck 
Observatory, in particular, Gregory Wirth. The analy- 
sis pipeline used to reduce the DEIMOS data was de- 
veloped at UC Berkeley with support from NSF grant 
AST-0071048. KSL thanks Eros Vanzella and Tucker 
Jones for sharing their data and for insightful comments. 
The authors thank Alice Shapley for sharing her com- 
posite spectra and for useful discussions. The spectrum 
of MS1512 cB58 presented here was kindly provided by 
Chuck Steidel. This work is based in part on observa- 
tions made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California 
Institute of Technology. We are grateful to the expert 
assistance of the staff of Kitt Peak National Observatory 
where the optical and near-infrared observations of the 
NDWFS Bootes Field were obtained. The authors thank 
NOAO for supporting the NOAO Deep Wide-Field Sur- 
vey. The authors also wish to recognize and acknowledge 
the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indige- 
nous Hawaiian community. 



REFERENCES 



Adelberger, K. L., Steidel, C. C, Pettini, M., Shapley, A. E., 
Reddy, N. A., Sz Erb, D. K. 2005, ApJ, 619, 697 

Adelberger, K. L., Steidel, C. O, Shapley, A. E., & Pettini, M. 
2003, ApJ, 584, 45 

Ashby, M. L. N., Stern, D., Brodwin, M., Griffith, R., Eisenhardt, 
P., Kozlowski, S., Kochanek, C. S., Bock, J. J., Borys, C, 
Brand, K., Brown, M. J. I., Cool, R., Cooray, A., Croft, S., 
Dey, A., Eisenstein, D., Gonzalez, A. H., Gorjian, V., Grogin, 
N. A., Ivison, R,. J., Jacob, J., Jannuzi, B. T., Mainzer, A., 
Moustakas, L. A., Rottgering, H. J. A., Seymour, N., Smith, 
H. A., Stanford, S. A., Stauffer, J. R., Sullivan, I., van Breugel, 
W., Willner, S. P., & Wright, E. L. 2009, ApJ, 701, 428 

Bouchet, P., Lequeux, J., Maurice, E., Prevot, L., & 
Prevot-Burnichon, M. L. 1985, A&A, 149, 330 

Bouwens, R. J., Illingworth, G. D., Blakeslee, J. P., Broadhurst, 
T. J., & Franx, M. 2004, ApJ, 611, LI 

Bouwens, R. J., Illingworth, G. D., Franx, M., Chary, R., Meurer, 
G. R., Conselice, C. J., Ford, H., Giavalisco, M., & van 
Dokkum, P. 2009, ApJ, 705, 936 

Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, 
ApJ, 670, 928 

Bouwens, R. J., Illingworth, G. D., Oesch, P. A., Franx, M., 
Labbe, I., Trenti, M., van Dokkum, P., Carollo, C. M., 
Gonzalez, V., Smit, R., & Magee, D. 2012, ApJ, 754, 83 

Calzctti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, 
J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 

Capak, P. L., Riechers, D., Scoville, N. Z., Carilli, C, Cox, P., 
Neri, R., Robertson, B., Salvato, M., Schinnerer, E., Yan, L., 
Wilson, G. W., Yun, M., Civano, F., Elvis, M., Karim, A., 
Mobasher, B., & Staguhn, J. G. 2011, Nature, 470, 233 

Chen, Y.-M., Tremonti, C. A., Heckman, T. M., Kauffmann, G., 
Weiner, B. J., Brinchmann, J., & Wang, J. 2010, AJ, 140, 445 

Cooper, M. C, Newman, J. A., Davis, M., Finkbeiner, D. P., & 
Gerke, B. F. 2012, Astrophysics Source Code Library, 3003 

Cooper, M. C, Newman, J. A., Weiner, B. J., Yan, R., Willmcr, 
C. N. A., Bundy, K., Coil, A. L., Conselice, C. J., Davis, M., 
Faber, S. M., Gerke, B. F., Guhathakurta, P., Koo, D. C, & 
Noeske, K. G. 2008, MNRAS, 383, 1058 

Daddi, E., Dannerbauer, H., Stern, D., Dickinson, M., Morrison, 
G., Elbaz, D., Giavalisco, M., Mancini, C, Pope, A., & 
Spinrad, H. 2009, ApJ, 694, 1517 

Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785 



Elbaz, D., Daddi, E., Le Borgne, D., Dickinson, M., Alexander, 

D. M., Chary, R., Starck, J., Brandt, W. N, Kitzbichler, M., 
MacDonald, E., Nonino, M., Popesso, P., Stern, D., & Vanzella, 

E. 2007, A&A, 468, 33 

Faber, S. M., Phillips, A. C, Kibrick, R. I., Alcott, B., Allen, 
S. L., Burrous, J., Cantrall, T., Clarke, D., Coil, A. L., Cowley, 
D. J., Davis, M., Deich, W. T. S., Dietsch, K., Gilmore, D. K., 
Harper, C. A., Hilyard, D. F., Lewis, J. P., McVeigh, M., 
Newman, J., Osborne, J., Schiavon, R., Stover, R. J., Tucker, 
D., Wallace, V., Wei, M., Wirth, G., & Wright, C. A. 2003, in 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 4841, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, ed. 
M. Iye & A. F. M. Moorwood, 1657-1669 

Ferguson, H. C, Dickinson, M., Giavalisco, M., Kretchmer, C, 
Ravindranath, S., Idzi, R., Taylor, E., Conselice, C. J., Fall, 
S. M., Gardner, J. P., Livio, M., Madau, P., Moustakas, L. A., 
Papovich, C. M., Somerville, R. S., Spinrad, H., & Stern, D. 
2004, ApJ, 600, L107 

Finkelstein, S. L., Papovich, C, Salmon, B., Finlator, K., 
Dickinson, M., Ferguson, H. C, Giavalisco, M., Koekemoer, 
A. M., Reddy, N. A., Bassett, R., Conselice, C. J., Dunlop, 
J. S., Faber, S. M., Grogin, N. A., Hathi, N. P., Kocevski, 
D. D., Lai, K., Lee, K.-S., McLure, R. J., Mobasher, B., k, 
Newman, J. A. 2012, ApJ, 756, 164 

Finkelstein, S. L., Rhoads, J. E., Malhotra, S., & Grogin, N. 2009, 
ApJ, 691, 465 

Gawiser, E., van Dokkum, P. G., Herrera, D., Maza, J., 

Castander, F. J., Infante, L., Lira, P., Quadri, R., Toner, R., 
Treister, E., Urry, C. M., Altmann, M., Assef, R., Christlein, 
D., Coppi, P. S., Duran, M. F., Franx, M., Galaz, G., Huerta, 
L., Liu, C, Lopez, S., Mendez, R., Moore, D. C, Rubio, M., 
Ruiz, M. T., Toft, S., & Yi, S. K. 2006, ApJS, 162, 1 

Giavalisco, M. 2002, ARA&A, 40, 579 

Giavalisco, M., Dickinson, M., Ferguson, H. C, Ravindranath, S., 
Kretchmer, C, Moustakas, L. A., Madau, P., Fall, S. M., 
Gardner, J. P., Livio, M., Papovich, C, Renzini, A., Spinrad, 
H., Stern, D., & Riess, A. 2004a, ApJ, 600, L103 



20 



Lee et al. 



Giavalisco, M., Ferguson, H. C, Koekemoer, A. M., Dickinson, 
M., Alexander, D. M., Bauer, F. E., Bergeron, J., Biagetti, C, 
Brandt, W. N., Casertano, S., Cesarsky, C, Chatzichristou, E., 
Consclicc, C, Cristiani, S., Da Costa, L., Dahlen, T., de Mcllo, 

D. , Eisenhardt, P., Erben, T., Fall, S. M., Fassnacht, C, 
Fosbury, R., Fruchtcr, A., Gardner, J. P., Grogin, N., Hook, 

R. N., Hornschcmcicr, A. E., Idzi, R., Jogee, S., Kretchmer, C, 
Laidler, V., Lee, K. S., Livio, M., Lucas, R., Madau, P., 
Mobasher, B., Moustakas, L. A., Nonino, M., Padovani, P., 
Papovich, C, Park, Y., Ravindranath, S., Renzini, A., 
Richardson, M., Riess, A., Rosati, P., Schirmer, M., Schreier, 

E. , Somerville, R. S., Spinrad, H., Stern, D., Stiavelli, M., 
Strolger, L., Urry, C. M., Vandame, B., Williams, R., & Wolf, 
C. 2004b, ApJ, 600, L93 

Gonzalez, V., Labbe, I., Bouwens, R. J., Illingworth, G., Franx, 

M., & Kriek, M. 2011, ApJ, 735, L34+ 
Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Armus, L. 

2000, ApJS, 129, 493 
Heckman, T. M., Robert, C, Leitherer, C, Garnett, D. R., & van 

dcr Rydt, F. 1998, ApJ, 503, 646 
Hildebrandt, H., Pielorz, J., Erben, T., Schneider, P., Eifler, T., 

Simon, P., & Dietrich, J. P. 2007, A&A, 462, 865 
Jannuzi, B. T. & Dey, A. 1999, in Astronomical Society of the 

Pacific Conference Series, 111 
Jones, T., Stark, D. P., & Ellis, R. S. 2012, ApJ, 751, 51 
Komatsu, E., Smith, K. M., Dunkley, J., Bennett, C. L., Gold, 

B. , Hinshaw, G., Jarosik, N., Larson, D., Nolta, M. R., Page, 
L., Spergel, D. N., Halpern, M., Hill, R. S., Kogut, A., Limon, 
M., Meyer, S. S., Odcgard, N., Tucker, G. S., Weiland, J. L., 
Wollack, E., & Wright, E. L. 2011, ApJS, 192, 18 

Kornei, K. A., Shapley, A. E., Erb, D. K., Steidel, C. C, Rcddy, 
N. A., Pettini, M., & Bogosavljevic, M. 2010, ApJ, 711, 693 

Kornei, K. A., Shapley, A. E., Martin, C. L., Coil, A. L., Lotz, 
J. M., Schiminovich, D., Bundy, K., & Noeske, K. G. 2012, 
ApJ, 758, 135 

Kuiper, E., Venemans, B. P., Hatch, N. A., Miley, G. K., & 

Rottgering, H. J. A. 2012, MNRAS, 425, 801 
Kulas, K. R., Shapley, A. E., Kollmeier, J. A., Zheng, Z., Steidel, 

C. C, & Hainline, K. N. 2012, ApJ, 745, 33 

Law, D. R., Steidel, C. C, Erb, D. K., Larkin, J. E., Pettini, M., 
Shapley, A. E., & Wright, S. A. 2007, ApJ, 669, 929 

Lee, K., Giavalisco, M., Conroy, C, Wechslcr, R. H., Ferguson, 
H. C, Somerville, R. S., Dickinson, M. E., & Urry, C. M. 2009, 
ApJ, 695, 368 

Lee, K.-S., Alberts, S., Atlee, D., Dey, A., Pope, A., Jannuzi, 

B. T., Reddy, N., & Brown, M. J. I. 2012a, ApJ, 758, L31 
Lee, K.-S., Dey, A., Reddy, N., Brown, M. J. I., Gonzalez, A. H., 

Jannuzi, B. T., Cooper, M. C, Fan, X., Bian, F., Glikman, E., 

Stern, D., Brodwin, M., & Cooray, A. 2011, ApJ, 733, 99 
Lee, K.-S., Ferguson, H. C, Wiklind, T., Dahlen, T., Dickinson, 

M. E., Giavalisco, M., Grogin, N., Papovich, C, Messias, H., 

Guo, Y., & Lin, L. 2012b, ApJ, 752, 66 
Lee, K.-S., Giavalisco, M., Gnedin, O. Y., Somerville, R. S., 

Ferguson, H. C, Dickinson, M., & Ouchi, M. 2006, ApJ, 642, 63 
Madau, P. 1995, ApJ, 441, 18 
Martin, C. L. 2005, ApJ, 621, 227 
— . 2006, ApJ, 647, 222 

Masters, D. & Capak, P. 2011, PASP, 123, 638 
Newman, J. A., Cooper, M. C, Davis, M., Faber, S. M., Coil, 
A. L., Guhathakurta, P., Koo, D. C, Phillips, A. C, Conroy, 

C. Dutton, A. A., Finkbciner, D. P., Gerke, B. F., Rosario, 

D. J., Weiner, B. J., Willmer, C. N. A., Yan, R., Harker, J. J., 
Kassin, S. A., Konidaris, N. P., Lai, K., Madgwick, D. S., 
Noeske, K. G., Wirth, G. D., Connolly, A. J., Kaiser, N., 
Kirby, E. N., Lemaux, B. C, Lin, L., Lotz, J. M., Luppino, 
G. A., Marinoni, C, Matthews, D. J., Metevier, A., & 
Schiavon, R. P. 2012, ArXiv c-prints 

Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713 



Ouchi, M., Shimasaku, K., Okamura, S., Furusawa, H., 

Kashikawa, N., Ota, K., Doi, M., Hamabe, M., Kimura, M., 
Komiyama, Y., Miyazaki, M., Miyazaki, S., Nakata, F., 
Sekiguchi, M., Yagi, M., & Yasuda, N. 2004a, ApJ, 611, 660 

— . 2004b, ApJ, 611, 685 

Ouchi, M. et al. 2005, ApJ, 620, LI 

Papovich, C, Finkelstein, S. L., Ferguson, H. C, Lotz, J. M., & 

Giavalisco, M. 2011, MNRAS, 412, 1123 
Pentericci, L., Grazian, A., Fontana, A., Castellano, M., 

Giallongo, E., Salimbeni, S., & Santini, P. 2009, A&A, 494, 553 
Pentericci, L., Grazian, A., Fontana, A., Salimbeni, S., Santini, P., 

de Santis, C, Gallozzi, S., & Giallongo, E. 2007, A&A, 471, 433 
Pettini, M., Rix, S. A., Steidel, C. C, Adelberger, K. L., Hunt, 

M. P., & Shapley, A. E. 2002, ApJ, 569, 742 
Pettini, M., Steidel, C. C, Adelberger, K. L., Dickinson, M., & 

Giavalisco, M. 2000, ApJ, 528, 96 
Quider, A. M., Pettini, M., Shapley, A. E., & Steidel, C. C. 2009, 

MNRAS, 398, 1263 
Reddy, N., Dickinson, M., Elbaz, D., Morrison, G., Giavalisco, 

M., Ivison, R., Papovich, C, Scott, D., Buat, V., Burgarclla, 

D., Charmandaris, V., Daddi, E., Magdis, G., Murphy, E., 

Altieri, B., Aussel, H., Dannerbauer, H., Dasyra, K., Hwang, 

H. S., Kartaltepe, J., Leiton, R., Magnelli, B., & Popesso, P. 

2012a, ApJ, 744, 154 
Reddy, N. A., Pettini, M., Steidel, C. C, Shapley, A. E., Erb, 

D. K., & Law, D. R. 2012b, ApJ, 754, 25 
Rcddy, N. A. & Steidel, C. C. 2009, ApJ, 692, 778 
Reddy, N. A., Steidel, C. C, Fadda, D., Yan, L., Pettini, M., 

Shapley, A. E., Erb, D. K., & Adelberger, K. L. 2006, ApJ, 

644, 792 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005, ApJS, 160, 115 
Schacrer, D. & Vacca, W. D. 1998, ApJ, 497, 618 
Schwartz, C. M. & Martin, C. L. 2004, ApJ, 610, 201 
Shapley, A. E. 2011, ARA&A, 49, 525 

Shapley, A. E., Steidel, C. C, Pettini, M., & Adelberger, K. L. 

2003, ApJ, 588, 65 
Stark, D. P., Ellis, R. S., Bunker, A., Bundy, K., Targett, T., 

Benson, A., & Lacy, M. 2009, ApJ, 697, 1493 
Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., & Bunker, A. 2010, 

MNRAS, 408, 1628 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., 

& Pettini, M. 1999, ApJ, 519, 1 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., Erb, D. K., 

Reddy, N. A., & Pettini, M. 2005, ApJ, 626, 44 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., Pettini, M., 

Dickinson, M., & Giavalisco, M. 2000, ApJ, 532, 170 
— . 2003, ApJ, 592, 728 

Steidel, C. C, Erb, D. K., Shapley, A. E., Pettini, M., Reddy, N., 
Bogosavljevic, M., Rudie, G. C, & Rakic, O. 2010, ApJ, 717, 
289 

Steidel, C. C, Pettini, M., & Adelberger, K. L. 2001, ApJ, 546, 
665 

Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 
D. K., Reddy, N. A., & Hunt, M. P. 2004, ApJ, 604, 534 

Vanzella, E., Cristiani, S., Dickinson, M., Kuntschner, H., 

Nonino, M., Rettura, A., Rosati, P., Vernet, J., Cesarsky, C, 
Ferguson, H. C, Fosbury, R. A. E., Giavalisco, M., Grazian, 
A., Haase, J., Moustakas, L. A., Popesso, P., Renzini, A., 
Stern, D., & The Goods Team. 2006, A&A, 454, 423 

Vanzella, E., Giavalisco, M., Dickinson, M., Cristiani, S., Nonino, 
M., Kuntschner, H., Popesso, P., Rosati, P., Renzini, A., Stern, 
D., Cesarsky, C, Ferguson, H. C, & Fosbury, R. A. E. 2009, 
ApJ, 695, 1163 

Verhamme, A., Schaerer, D., & Maselli, A. 2006, A&A, 460, 397 
Weiner, B. J., Coil, A. L., Prochaska, J. X., Newman, J. A., 
Cooper, M. C, Bundy, K., Conselice, C. J., Dutton, A. A., 
Faber, S. M., Koo, D. C, Lotz, J. M., Rieke, G. H., & Rubin, 
K. H. R. 2009, ApJ, 692, 187 



